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Synopsis 


The machining conditions for a particular operation oyLStill to a large extent obtained from 
manufacturer’s data or machining handbooks which provide rather conservative estimates. 
These machining conditions are normally obtained on the basis of economic criteria such as 
minimum production cost or time per piece, or maximum profit rate. More realistic condi- 
tions can be obtained by considering economic optimum under technical constraints of real 
life machining such as maximum allowable cutting power, cutting force, surface roughness 
etc. Machining conditions thus obtained will be the optimal conditions. 

Because of non-availability of nonlinear programming (NLP) algorit hms , earlier efforts 
for finding optimal machining conditions have been limited to differential calculus based 
methods or numerical methods. The number of constraints considered were only a few 
because inclusion of large number of constraints involve development of complex logic for 
searching the optimum. Later NLP methods, both unconstrained and constrained, have 
been used for evaluating the optimal machining conditions. Geometric programming (GP) 
has also been extensively used for obtaining optimal machining condition since it was ob- 
served that a large number of technical constraints could be represented by monomials or 
single term polynomials. 

The real fruits of optimization of machining conditions have not been reaped because of 
two difficulties. (1) The unknown number of passes at multipass problem formulation stage 
and (2) the requirement of a feasible or non-feasible starting solution during optimization 
stage which again requires solution of another optimization problem. The first problem 
has been solved to some extent by taking a small depth of cut in finish pass and trying all 
possible number of even passes. Another solution is to solve the multipass turning problem 
using dynamic programming (DP) approach where each pass is considered a stage. In this 
case the total depth is divided into a number of segments and single pass problems are 
recursively solved to reach the optimal solution. The optimal solution contains the optimal 
number of passes alongwith the number of segments falling into each pass, and cutting 
speeds and feeds. Any effort to reduce the segment size squarely increases the number of 
single pass solutions. This method, therefore, requires more computation time and is not 



popular. For the solution of second problem sequential quadratic programming method 
(SQP) may be used. This method can be started from any arbitrary point but it is not 
always easy to implement. 

The present work has been carried out mainly for removing the above mentioned diffi- 
culties. The effort has been to find the optimal number of passes uniquely. The minimum 
production cost per piece criterion has been adopted for the determination of optimal ma- 
chining conditions. The multipass problem has not been solved by taking all passes together, 
instead the solution has been synthesized using single pass solutions. While attempting to 
get the solution the following points have also been considered:, 

• Selection of computationally efficient and easy to implement methodology for the solu- 
tion of single pass turning problem. The methodology selected should not necessarily 
require a feasible or a non-feasible starting point and should not depend upon the 
optimization parameters. 

• Parametric analysis of single pass methodology so that solution to new problems 
generated because of changed workpiece dimensions can be easily computed. This 
will help in reduction of computation time during application of DP approach to 
multipass turning problem. 

• A method for logical selection of depth of cut for finish pass and optimal distribution 
of depth of cut during finish as well as rough passes. 

• The concept of uneven passes as much of the research work is still concentrated towards 
even passes. 

• Development of a model for optimal clubbing of machinable volumes, horizontally, 
vertically or in circular order for achieving minimum production cost. 

A detailed review of existing literature is presented in Chapter 1. The topics covered in- 
clude CAPP systems for rotational parts, various tool life equations, different optimization 
criteria, optimization techniques, etc. The papers have also been summarized in a tabular 
form which indicates important factors such as single or multipass turning, economic crite- 
rion, deterministic or probabilistic model, unconstrained or constrained turning constraints 
used in the model, optimization approach, solution methodology, decision variables and 
type of tool life equation used for constructing the tmning model. 

The various time components viz. preparation time, idling time, machining time, tool 
changing time and the corresponding costs have been defined in Chapter 2. The technical 
constraints such as cutting force, cutting power, surface finish, etc. have also been discussed 
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in this chapter. It has been observed that majority of constraints are monomials. There- 
fore, a general multipass cost minimization problem has been formulated using extended 
Taylor’s tool life equation and monomial constraints. Assuming that one finish and a few 
even rough passes are required to remove total material, the minimum and maximum total 
number of passes are decided on the basis of allowable ranges of depths of cut for finish 
as well as rough passes. The multipass problem is repeatedly solved for all possible total 
numbers of passes and the niunber giving the least minimum production cost or time per 
piece is declared optimal. The solution process is further simplified by providing depth of 
cut distribution also along with number of passes as input to the multipass problem. The 
common practice is to use a small depth of cut for finish pass and remove the remaining 
depth in even member of rough passes. The multipass problem is now decomposed into a 
number of single pass problems which are easier to solve. At this stage assumptions of equal 
workpiece diameter and/or equal tool fife in all passes are also used to reduce the number of 
single pass problems. The general multipass turning^may also be solved using the dynamic 
programming (DP) approach which gives urdquely optimal number of passes using single 
pass solutions. The total depth of cut is divided into a number of segments and the optimal 
solution contains optimal number of passes, the number of segments falling into each pass, 
and the cutting speed as well as feed for each pass. 

Out of two approaches, optimizing all the passes simultaneously or using single pass so- 
lutions to synthesize results for multipass turning problem, the later one has been adopted 
in the present work as it is easier to obtain single pass solutions. The common practice is 
to search the optima in speed-feed plane. Various NLP methods use different techniques 
to formulate and update search vector comprising of speed and feed components and then 
use single dimensional search methods to obtain optima along these search vectors. It has 
been found that at a given tool life single pass production cost for a given depth of cut 
is minimized by using maximmn allowable feed under constraints at that particular tool 
life. Therefore, instead of finding a search vector and updating it in speed-feed plane, a 
one dimensional search along tool life may be performed to find out the optimal tool life 
and feed corresponding to minimum production cost. The cutting speed is obtained using 
extended Taylor’s tool fife equation. Thus, the process of searching for optima under the 
constraints is immensely simplified. 

The single pass solution methodology has been described in Chapter 3. Parametric 
analysis for this methodology has also been discussed. This analysis helps in finding the 
solution for new single pass truning problems with changed workpiece dimensions for the 
same depth of cut without solAung the new problem. 
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In this chapter two solution models have been presented. The first one uses a rather 
new concept of assumed maximum number of passes and depth of cut verses minimtun 
production cost trends. The minimmn production costs for a series of feasible depths of cut 
are obtained for finish as well as rough passes. One finish pass and a large niunber of rough 
passes are assumed to remove the total depth of cut. A binary V 9 .riable is associated with 
each feasible depth in series in finish as well as rough passes. A linear integer programming 
model is formed using the feasible depths in series, the minimiun production costs corre- 
' sponding to the depths, and binary variable for the assumed maximiun number of passes. 
The problem is solved using LINDO optimizer and optimal number of passes axe obtained 
by counting non-zero binary variables in the solution. 

The second model is based on dynamic programming approach and is already available 
in the literature. It could not be explored fully because any effort to reduce the segment size 
squarely increases the number of single pass solutions. The parametric analysis for single- 
pass solution methodology suggested in this chapter reduces the computational work and 
thus helps in reduction of segment size to 0.1 mm. Earlier workers had taken the segment 
size as 1 mm. 

The numerical results for a few straight multipass tinning examples have been presented 
in Chapter 4. It has been observed during the use of integer programming model that un- 
even passes yield lower minimum production cost. DP approach gives optimal solution 
with better accuracy and is capable of handling the general multipass problem without 
making assumptions such as equal workpiece diameter or equal tool life in all passes. DP 
approach along with single, pass solution methodology has been comipared with Sequential 
Unconstrained Minimization Technique (SUMT) with Davidon Fletcher Powell’s method 
(DFP), Generalized Reduced Gradient (GRG) method and Sequential Quadratic Program- 
ming (SQP) method. It has been observed that for a general multipass turning example, DP 
approach is easier to implement than any other method and it does not require a feasible or 
nonfeasible starting point. Further, this approach requires less computer time than SUMT 
with DFP and GRG and works equally well as SQP. 

\ 

Rough turning of stepped shaft involves decision on horizontal, vertical, or circular club- 
bing of machinable volumes. A model has been suggested to perform logical clubbing of 
volumes in Chapter 5. This model has been solved using LINDO optimizer. The results 
obtained for two shaft which are similar in shape but of different dimensions show that 
different optimal grouping may be obtained even if the shape of shafts appear to be similar. 
Further, the optimal clubbing of voliunes obtained using mathematical model yields better 
solution than that obtained using horizontal clubbing, or using a heuristic of ascending 
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mininmm production cost/unit volume. 

The conclusions and scope for future work have been presented in Chapter 6. There 
are two basic issues involved in the solution of multipass turning problems. The first is to 
obtain optimal number of passes and second is to treat all the passes simultaneously, or 
decompose the multipass problem into single pass problems. It is easier to solve single pass 
problems. Therefore, second issue has been taken care by suggesting a simple single pass 
solution methodology. The parametric analysis for this methodology has been performed 
for the first time. A new integer progra mming model has been suggested for obtaining the 
optimal number of passes as well as the optimal depth of cut distribution. The available 
DP model has been adopted because it uniquely gives optimal number of passes. During 
application of dynamic programming approach to multipass turning operation the current 
practice is to optimize the single pass turning problems each time they are formed due to 
varying workpiece diameter even if the depth of cut remainsthe same. Once solution to a 
single pass turning problem for a given depth of cut is available the parametric analysis 
helps in getting the solution of these newly formed problems of same depth of cut with 
varying workpiece diameter in each pass through simple calculations. Both the integer pro- 
gramming model and DP approach have been useful in logically deciding the depth of cut 
for finish pass which used to be kept at some small or least value during the finish pass. 

' Further, these models have helped in obtaining better minimum production cost by explor- 
ing uneven passes. Parametric analysis for single pass solution method has been very useful 
in reducing the segment size of depth of cut and computational effort during DP approach. 

The present work considers monomial constraints and may be extended to explore the 
treatment of nonmonomial constraints. The turning operation has been taken for simplicity. 
It will be interesting to know if other machining operations such as milling, drilling, reaming, 
boring, etc. can also be treated using the same methodology. 
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Chapter 1 

Introduction and Literature Review 


1.1 INTRODUCTION 

The world of machining has been changing continuously and rapidly. Most of the earlier 
changes were brought about by advancements in machining technology and development 
of new tool materials. Rapid changes in the past few decades can be attributed mainly to 
automation leading to the development of programmable machining centers like NC, CNC, 
DNC, etc. These equipped with supporting activities such as tool changing, automatic work 
loading/unloading, onUne process monitoring and control, etc. have significantly enhanced 
machining capabilities. Such developments in machining have been motivated more by the 
need to remain competitive in global market. The major effect of these changes can be seen 
in the increased productivity through substantial reduction of idle time for machines as well 
as parts. Because of automation and monitoring, the effective operating time - uptime for 
machines has gone up from 10 % to 65 % (Fig. 1.1). With the improvement in operating 
time, much more rewards can be reaped from the optimization of machining conditions 
than with earlier machining environment where uptime itself was quite low. Yet another 
factor necessitating the optimization of machining conditions is capital intensive investment 
required by automated machines. Figure 1.2 clearly shows that cost vs cutting speed curves 
are much steeper for CNC machines than those for conventional machines. Consequently, 
the optimal cost comes down with the use of NC/CNC machines but the cost of being away 
from the optimal cutting speed increases rapidly (Fig. 1.2b). Hence it becomes imperative 
to work at optimal machining conditions as far as possible. 

Technical limitations can be critical during machining operations. These limitations 
usually relate to lack of requisite power and torque, excessive force, bad surface finish, and 
limitations on speed, feed and depth of cut. As proposed by Lundholm [1991], the machining 
conditions may be classified as technical, recommended, global economical optimal, and 
optimal. Their interrelationship is displayed on feed - cutting speed plane in Figure 1.3. 
The feasible region defined by machining conditions (speed, feed, depth of cut and total 
number of passes) is bounded by various constraints. A feasible region therefore contains 
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Figure 1.1: Economic importance of monitoring in manufacturing process [Tonshoff et al., 
1988]. 
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Figure 1.2: Production costs per piece as a function of the cutting speed: (a) conventional 
machining, (b) CNC machining [Tonshoff et al. 1988]. 



all possible values of machining conditions satisfying a given set of constraints. 

Technical Machining Conditions (TMC) 

As shown in Figure 1.3, TMC encompass the whole set of feasible values of machining 
conditions satisfying all the given constraints in the form of limitations on cutting force 
and cutting power, surface finish requirements, practical bounds on speed, feed and depth 
of cut, geometrical and dimensional parameters of the workpiece, and tool and workpiece 
materials. 



Figure 1.3: Classification of machining conditions [Lundholm, 1991]. 

Recommended Machining Conditions (RMC) 

These are obtained in the form of recommendations from data handbooks [Metcut, 
1980; Ai, 1985] which have been compiled over the years on the recommendations of tool 
suppliers, the experience of manufacturers (users), synthetic analysis of machining processes, 
and workers’ experience. RMC values are invariably conservative. Various expert systems 
have been developed based on RMC to provide recommended machining conditions. One 
such system provides the recommendations which are stored in a computer and are retrieved 
during automated process planning [Wang and Wysk, 1986]. 

Global Economical Optimal Machining Conditions (GEOMC) 

These machining conditions are calculated taking only the economical aspect of the pro- 
cess into account, while completely disregarding the technical limitations. In other words. 



GEOMC are determined using unconstrained optimization with the objective comprising 
various manufacturing costs. Therefore, under many situations, GEOMC lead to unreason- 
ably high values of speed, feed and depth of cut. 

Optimal Machining Conditions (OMC) 

These are obtained considering economic optimum under technical limitations. OMC, 
in general, are more favourable than TMC, RMC and GEOMC. These can be determined 
online as well as offline. 

The online determination of OMC requires reliable initial setting of speed, feed and 
depth of cut. These input machining conditions are optimized online during machining of 
the workpiece. The online OMC determination involves sensing and monitoring of various 
events viz. excessive motor current, worn tool, chipping, tool breakage, cutting forces, 
catastrophic chatter, etc. In order to sense and monitor each of these events separately, or 
in combination of other events, a niunber of computer based subsystems are developed. All 
of these subsystems transmit necessary information to a main controller which is governed by 
a master computer. An optimization algorithm determines the online OMC after sensing the 
real time data from various subsystems. The optimal machining conditions thus obtained 
are supplied to the machining center. This cycle of sensing and monitoring, and optimization 
of machining conditions is repeated in real-time. 

Often all machining conditions cannot be optimized on line. For example in turning 
operations, the cutting speed and depth of cut are optimized offline and are kept constant 
during actual machining, whereas the feed is optimized online to obtain OMC. A more 
useful system shall be that which optimizes both speed as well as feed online keeping depth 
of cut as constant [Lundholm, 1991]. 

The offline determination of OMC is performed dining planning stage. It may be in- 
tegrated with machine selection, tool material and geometry selection, and operation se- 
quencing. Some of the recent computer assisted process planning (CAPP) systems have 
started using mathematical optimization for offline OMC determination. Table 1.1 sum- 
marises some of the popular CAPP systems using OMC determined offline. Offline OMC 
may also act as reliable input to onUne OMC determination. 

1.1.1 Tool Selection 

The tool selection essentially means the selection of (i) tool holder and (ii) geometry and 
material of insert. The selection of tool holder depends on the tool post on one hand and 
the size and type of insert on the other. Insert selection depends upon the geometrical and 
dimensional features of blank and final finished part as well as on the machining condition 
i.e., ranges of speed, feed and depth of cut. The machining conditions, on the other hand, 
can be determined only when the features of insert are known. This difflculty is, however, 


overcome by the fact that tool selection, does not require the exact value of machining data 
but only their ranges [Giusti et al., 1986]. Thus, after selecting the tool the exact machining 
conditions may be computed considering the economic factors expected to prevail during 
machining. 

Considering the above the tool selection becomes a two stage process. In the first stage a 
list of tools is generated on the basis of knowledge of cutting technology. It takes into account 
the problems of geometrical compatibility with workpiece, and differences between roughing 
and finishing tools. The choice of best tool holder and the insert is not an easy task. It is 
because the choice criterion often can not be expressed as an exact knowledge. This task gets 
further complicated because of continuous additions of new tool materials and geometry. 
An automatic tool selection system should therefore be such that it can be updated without 
a software specialist. An expert system based on algorithmic if (conditions)-then (action 
solutions) rules can be a good solution to the problem of tool selection [Guisti et al., 1986]. 
The other solution method for the first stage is called production rule matrix method and is 
based on production rules represented in a tabular foxm instead of logical (if-then) statement 
form [Domazet, 1990]. This method uses integer multiplication to select a tool and therefore 
works faster than the expert system. The tool selected at this stage may not be optimal and 
therefore the second stage is used to find the optimal tool from the list of tools generated 
at the first stage. 

In the second stage, an economic criterion is used to find OMC corresponding to a tool 
from the list of tools selected in the first stage. The optimal tool from the list of tools may 
be found using an exhaustive search over a number of tools using the selected criterion. 
Chen et al. [1989] have presented a heuristic method to automatically select cutting tools 
for rough turning operations. The selection of tool from the appropriate list of tools is 
based on an economic objective viz. minimum production cost. The cost of machining with 
a given tool is estimated after the determination of cutting conditions imder constraints. 
It is possible to estimate, after examination of constraints, whether the next tool in the 
list will improve cutting conditions and will give reduced minimum production cost or not. 
The heuristic method works faster and eliminates the need for exhaustive search. The tool 
selected using the heuristic method may not always be optimal. In case of any difficulty in 
the implementation of first stage, a list of already available tools in the library may also be 
used as input to the second stage. 

In case only tool material is to be selected (assuming tool geometry to be known), any 
method to determine OMC can be used and the tool material corresponding to the best 
value of economic criterion is selected. Lukic et al. [1991] present the Yugoslav approach 
to solve this problem. This approach makes use of constants of Taylor’s tool life equation 
and is based on some elaborate experimentation. Fenton and Gagon [1993] have plotted 
the minimum production cost per piece (Qp) and minimum production time per piece {tfo) 


K 



points on two axis of a graph for different tool materials. The combination of minimum 
production cost and minimum production time points which is nearer to the origin is selected 
and the tool is required to work between these two points. 

1.1.2 Tool Life 

In a pioneering work, Taylor [1907] suggested a relationship between machining time and 
tool life for achieving minimum production cost. The relationship has come to be known 
as Taylor’s tool life equation and relates cutting speed (v) and to tool life (T) as 

vT^^=Cto ( 1 - 1 ) 

where ni and Cto are constants. 

The above equation proposed nearly a century ago is still widely used [Lindstrom, 1989]. 
The relationship is valid for a given value of feed and depth of cut. An extended form of 
the above equation which includes the effect of feed (/) and depth of cut (d) is given as 
[Armerago and Brown, 1969], 

Exponents ni, n 2 , and constant Ct niay be evaluated from test data using response surface 
methodology, ARMA modelling, regression analysis or neural network [Chryssolouris et al., 
1990]. The number of tests may be significantly reduced by using random strategy method 
[Kuljanic, 1980]. Another form of the above equation defining tool life is 

y Ko 

v“o/A)(i'ro 

where 

Ko = ; ao = (l/ni) ; / 3 o = (ua/ni) ; 70 = (n^/rii) 

Golding [1959] suggested a tool life equation in terms of equivalent chip thickness (he) 
and is of the form, 

= Qri (1.3) 

where Cti is a constant and exponent rii is a function of v and hg and exponent is a 
function of only. Evaluation of these exponents is not simple, therefore eqn. (1.3) is 
usually not used for OMC determination. Brewer and Reuda [1963], however, used this 
equation for OMC determination treating and as constants for a given tool-workpiece 
combination over some working range of v and he- Based on the results of machining eight 
different types of cast irons with three different t 3 rpes of tool materials (HSS, carbide and 
ceranaic), they concluded that exponent 714 generally depends only on tool material while 
exponent Uq depends only on work material. In other words, exponents n 4 and nj are 
constants for a given tool-work material and are independent of v and he- The equivalent 
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chip thickness [he] (a ratio of area of cut to contact lenth is defined as [Brewer and 
Reuda, 1963] 


Hence eqn. (1.3) becomes 


K = Gf 


vT^^iGf)^ = Cti 


(1.4) 

(1.5) 


or 


and to a very close accuracy 


G 




d 

d-r{l-sini>) , ( 90 — 1 />) 

COSTj) ‘ 180 


( 1 . 6 ) 


(1.7) 


Here r and represent the nose radius and side cutting-edge angle of tool, respectively 
(Fig. 1.4). 



Figure 1.4: Chip equivalent concept [Brewer and Reuda, 1963]. 

Yellowley and Gunn [1989] have approximated G as 

G = {d) / (n^d n-j) (f-8) 

where and n-^ are constants. Approximating G = d"® in eqn. (1.8) leads to the extended 
tool life eqn. (1.2 ). From eqns. (1.2) and (1.6), it can be easily seen that ni = n^\ na = ns; 
^3 = (nsns) and Ct = Cti- 
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Exponents ni, n 2 , may have elements of uncertainty and hence, tool life may turn 
stochastic. However, it is found that the deterministic models are fairly close to their 
stochastic cotmterparts in prescribing OMC. Since eqn. (1.3) has been approximated as 
eqn. (1.2) and it is widely used for OMC determination. The tool' life given by eqn. (1.2 ) 
has been used in the present research work. Exponents rij, n 2 and have been assumed 
to be constant, and deterministic tool life equation has been used to construct production 
cost and production time functions. 

1.1.3 Technical Constraints 

In real life, turning operation is restricted by various practical constraints in the form of 
cutting power, cutting force, spindle torque, minimum and maximum tool life, tool thrust, 
tool wear, maximum tensile stress on tool rake face, maximiun temperature in primary shear 
zone, temperature at chip-tool interface, continuous chip flow requirement, chip breaking 
region requirement, surface finish and surface integrity requirements, allowable limits of 
speed, feed, depth of cut, diametral accmacy and any other related constraint. These con- 
straints are represented as functions of machining conditions elements speed, feed and depth 
of cut. Their functional relationships shall be presented and used in OMC determination 
in section 2.3 of the next chapter. 

1.1.4 Economic Criteria 

Machining conditions affect production time and tool life and hence, the production cost. 
The problem of OMC determination is formulated on the basis of some such criterion. The 
following two criteria are most commonly used ; 

1. Minimum production cost: OMC are decided such that cost of producing a workpiece is 
minirmim - 

2. Maximum production rate or minimum production time: If the operation is a bottleneck 
in production line, it is required to be produced in the least possible time. 

Apart from these, other criteria in vogue are: (a) maximxun tool life, (b) maximmn profit 
rate, (c) minimum of weighted average of production cost and production time, (d) maxi- 
mum material removal rate and (e) maximum material removal per unit flank wear. 

The list is not exhaustive. The criterion of production time per piece, is equivalent to 
the criterion of production cost per piece. Therefore the method which can be applied to 
minimization of production cost can also be applied to minimization of production time. 
The tool life requirement, either on time scale or on the basis of number of pieces produced, 
may be treated as a constraint during minimization of production cost or production time. 
Therefore, in the present work, the minimum production cost criterion has been adopted 
out of the criteria listed above. 
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1.1.5 Optimization Procedure 

Optimization is the act of obtaining best results under given circumstances defined by cri- 
terion and constraints. As the effort required and benefit accrued in any practical situation 
may be expressed as a function of certain variables, optimization, is the process which gives 
maximum or minimmn of one or more objective functions subject to a set of constraints 
and thus guides in the choice of values of decision variables. 

Decision Variables 

Any optimization problem is conceived with the help of interrelated quantities out of which 
some are viewed as variables during the design process. In general, certain quantities are 
fixed at the outset and are called fixed input parameters and some others are treated as 
decision variables. In the problem of OMC determination for multipass turning, speed, feed, 
and depth of cut alongwith ntunber of passes, in general, form the set of decision variables 
and various machining constants and mechanical and metallurgical properties are treated as 
constants. Sometimes, in order to make the model simple, a few of these decision variables 
are fixed during the input stage to make them fixed parameters. 

Constraints 

In many problems, decision variables caimot be chosen arbitrarily; rather they have to satisfy 
certain specified functional and other requirements. The restrictions which must be satisfied 
in order to produce an acceptable decision are collectively called decision constraints. The 
technical constraints described in section 1.1.3 are the constraints in the problem of OMC 
determination. 

Objective Function 

Objective function is a measure which guides in the choice of optimal design out of a 
number of acceptable decisions available. It is a criterion used for comparison of different 
acceptable decisions for selecting the optimal one. This criterion, when expressed as a 
mathematically computable function of decision variables is known as objective function. 
The choice of this function is governed by the nature of problem. The various economic 
criteria for optimization purpose have been described in section 1.1.4. An objective function 
can be formulated on the basis of any given criterion. As stated ear her, the criterion used 
for optimization in the present work is Tninimum production cost because minimization of 
production time is similar to minimization of production cost. The tool life requirement 
has been included as a constraint. 





1.2 Literature Survey 


The important factors over which the research reported may be distributed are: (a) eco- 
nomic criterion, (b) single or multipass machining, (c) type of model (probabilistic or de- 
terministic), (d) the solution space (unconstrained or constrained; if constrained then type 
of constraints), (e) optimization approach, (f) solution methodology (use of different as- 
sumptions to simplify the model), (g) fixed parameters (d, / and/or n), (h) optimal output 
decision variables and (i) t 3 ^e of tool life equation (extended Taylor’s or modified Gold- 
ing’s). A su mm ary of available literature in the field of optimization of single or multipass 
turning is presented in Table 1.2. Five different t 3 rpes of models have been used for modelling 
turning of a cylindrical workpiece. These are discussed in subsequent subsections. 


1.2.1 Single Pass Turning Models 


This is the basic model for turning. It assmnes a constant workpiece diameter and a constant 
depth of cut. The production cost (Cto) or production time (tto) fimctions when partially 
differentiated with respect to speed or feed and equated to zero do not give unique feed and 
speed. According to Armerago and Brown [1969], if either feed or speed is fixed, a point 
corresponding to minimum production cost or production time can be obtained. There exist 
a tool life (Tmc) for a given depth of cut and feed corresponding to nainimum production 
cost and is obtained as 

= (ao - (1.9) 

Ki 

and the speed (v^c) corresponding to the above tool life T^c is 


Kpi 


( 1 . 10 ) 


where ki, kt and tc are over head cost per unit time, tool cost per edge and tool changing 
time per edge respectively. Similarly tool life (Tmt) corresponding to minimum production 
cost is given by 

Trat = (ao - l)tc (1-11) 


and the speed (v^t) corresponding to the above tool life is 




Kq^ 

{Trntf^d-ro)^ 


( 1 . 12 ) 


The earlier efforts (Table 1.2) for OMC determination have been mostly limited to the 
selection of speed using these equations and for a given feed fixed on the basis of experience. 
In constrained situations if the resulting speed violates any of the constraint then speed is 
lowered till that constraint is satisfied. If the constraint is not satisfied even at the lowest 
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possible speed then the next lower feed is selected and again the optimum cutting speed is 
found using eqn. (1.10) or (1.12) [Brewer and Reuda, 1963; Field et al., 1968]. 

The research work started with the consideration of all possible combinations of speed 
and feed and selection of the optimum cutting speed and feed corresponding to minimnni 
production cost or minimum production time for unconstrained turning. Crookal [1969] 
has solved the single pass turning problem with constraints using performance envelope 
concept. It is a graphical approach and many graphs have to be plotted in order to find 
the cutting speed and feed through performance envelope. This method although could not 
become popular, a good concept viz. production cost or production time per unit turned 
length has been used for the first time in this work. 

The constrained single pass turning model has been solved using different nonlinear 
programming (NLP) methods viz. Lagrange Function Method [Bhattacharya et al., 1970], 
Sequential Unconstrained Minimization Technique (SUMT) using Newton Raphson method 
or Davidon Fletcher and Powell (DFP) method [Iwata et al., 1972; Hati and Rao, 1976], 
Classical Gradient Method [Elkarmany and Papal, 1978], Neldermead Simplex Method 
[Agapiou, 1992a & b]. The last method has been used for weighted average of production 
cost and production time functions. 

Since all constraints except the constraints of tool wear (section 2.3.8) and chip-breaking 
region (section 2.3.12) are monomials, the geometric programming (GP) has been exten- 
sively used for minimization of production cost or production time functions derived using 
extended tool life eqn. (1.2) for single pass turning with constraints. Cto or tto function with 
one constraint constitutes zero degree of difficulty problem (three terms and two variables). 
Any additional constraint introduced to the problem increases the degree of difficulty by 
one. Single degree difficulty problem has been solved through maximization of dual prob- 
lem using differential calculus [Walvekar and Lambert, 1970]. Petropoulus [1973] has solved 
this problem through numerical iteration. Ermer and Kromodihardjo [1981] have solved the 
problem of minimum production cost or production time with any number of constraints 
(a general degree diffictilty problem) using GP-LP method [Kochenberger et al. 1973]. The 
minimization of a posynomial function under posynomial constraints can be performed us- 
ing this method. However, the rate of convergence for the method becomes slow if more 
constraints are added [Averiel, 1980]. Eskcioglu [1985] has attempted to solve this Geomet- 
ric Programming Problem (GPP) by analysis. The two variable problem is solved for all 
combinations of zero and one degree difficulty problems. The maximum of the dual is se- 
lected from all zero and one degree difficulty problems and corresponding primal variables, 
i.e., speed and feed are found. If this speed-feed combination lies in^the feasible range of 
speed and feed, it is declared as the solution. In fact, it is not the correct solution because 
the solution which lies in the feasible range of only speed and feed may still lie outside the 
feasible region defined by all the constraints. Recently Gopalkrishnan and Faiz-AlKhyyal 
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[1991] have solved one degree difficulty problem by an analytical method. Two zero degree 
difficulty problems are solved uniquely in dual space and one single degree problem is solved 
taking intersection point of the two constraints in In v - In f space. The minimum point 
is located out of the three points using some extra logic. The development of extra logic is 
not an easy task when there are more than two constraints in single pass turning problems. 
However, this logic can be developed for such problems because the minima of the primal 
problem lies either on one constraint or on extreme point of the feasible region defined by 
intersection of all constraints. The feasible region for monomial constraints is a polygon in 
In v-ln f space. 

Okushima and Hitomi introduced the maximimi profit rate concept in 1964 [Hitomi 
1971]. Wu and Ermer [1966] presented a method to determine the maximum profit rate 
cutting speed for single pass turning case without constraints using the principle of marginal 
revenue and marginal cost. This model gives the value of optimiun speed for a given feed. 
Armerago and Russel [1966] have analysed the profit rate function for single pass tmning, 
treating speed and feed as variables. Through (dPr/dv) = 0 and {dPr/df) = 0 nomographs 
they have shown that no unique speed-feed point exists which satisfies these equations. 
They have also presented the nomographs for different constraints to show how and when 
these become active. Their attempts have been of theoretical nature only. Hitomi [1971] has 
presented a computational method to calculate speed while keeping the feed constant. The 
method has been developed on the basis of differential calculus and the unique combination 
of speed and feed is obtained using a non-linear programming (NLP) technique. Boothroyd 
and Rusek [1976] have presented the equation to calculate the tool life (Tpr) corresponding 
to maximum profit rate for a given revenue S and fixed feed and depth of cut as 



1-ni ktLIp 

ni '' ktS n^SvrTr 


(1.13) 


where for a reference speed and tool life , the tool life equation is given as 


VrTr^^ = CtO 


(1.14) 


The speed corresponding to the above tool life {Tpr) is obtained as 


V, 


'pr 


Ctq 

n^ni 

■^pr 


(1.15) 


These equations can be solved numerically and the tool life corresponding to maximum 
profit rate can be obtained for a given feed. Once tool life is known, the speed correspond- 
ing to this tool life can also be obtained. Hati and Rao [1976] have used sequential un- 
constrained minimization technique (SUMT) alongwith Davidon Fletcher Powell’s method 
(DFP) to maximize profit rate under constraints for single pass turning. Any general non- 
linear programming method may be applied to this kind of problem. However, GP which 
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found extensive use in minimization of production cost and production time could not be 
applied to the maximization of profit rate problem. The reason is that profit rate function 
is not posynomial [Ermer and Kromodihardjo, 1981]. Gupta et al., [1994] have converted 
the problem of maximization of profit rate function under posynomial constraints to min- 
imization of a posynomial function under posynomial constraints. It has been achieved 
by adding few variables and one posynomial constraint to the original problem and the 
resulting minimization problem is amenable to GP-LP technique. 

The knowledge of selling price S is essential for the determination of profit rate and the 
solution is sensitive to the value of S. Further, it may not be possible to determine the 
allocation of S for each operation on a workpiece and, therefore, the maximization of profit 
rate criterion is not commonly used as an economic criterion. 


1.2.2 Multipass Turning Models 

Multiple passes are adopted in cases where it is not possible to obtain the desired surface 
finish in a single pass, or it is not possible to remove the total stock, dji, in a single pass. 
For multipass turning, the number of passes Up and the depth of cut for individual passes 
di] i = 1, . . . , Up are important decision variables in addition to the speed and feed for 
individual passes. 

The total number of passes Up includes one finish pass with smaller depth of cut dg and n 
= (rip — 1) number of rough passes such that n^m < i^min and i^maxj fb® minimum 

and maximum possible nrunber of rough passes, can be evaluated as 


'^min 


dti dg 


(1.16) 


and 

(1.17) 

where drmm and drmax are minimum and maximum values for depth of cut in a rough 
pass, and [ ]'^ or [ ]“ imply integer on higher or lower side of the real value of term 
in brackets. Therefore, the diameter of the workpiece at the start of a pass is not the 
same for all passes as in multipass tmning operation. The models have been generally 
developed on the basis constant initial or final workpiece diameters ignoring the diameter 
changes during intermediate passes. The production cost (or production time) per piece 
in multipass turning is the sum of production costs for individual passes. The production 
cost based on the initial diameter will therefore be an upper botmd on minimum production 
cost, whereas production cost based on the final diameter will be a lower bound. The same 
is true for minimum production time function also. Better modelling, therefore, requires 
consideration of diameter changes dxiring each pass. 


dfi dg 
drmin 
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Research work in the area of multipass turning may be classified on the basis of selection 
of depths of cut for all passes. Important selection procedures are: (a) a small depth in the 
last pass and equal depths in rough passes, (b) unequal depths of cut decided in advance, (c) 
variable depths of cut and dynamic programming approach, (d) variable number of passes 
and depths of cut and (e) largest cut and heuristic for turning of stepped shaft. 

(a) A Small Depth in the Last Pass and Equal Depths in Rough Passes 

A general rip pass turning problem reduces to rip single pass turning problems if "small 
depth in the last pass (to get the desired surface finish) and equal depths in rough passes" 
strategy is adopted. Further, if constant workpiece diameter assumption is used in place 
of reduced workpiece diameter after each pass then all rough passes are the same and 
hence only two (one finish and one rough pass) single pass turning problems have to be 
solved. Earlier researchers used this approach to solve multipass turning problems. Berra 
and Barrash [1968] used empirical equations to obtain machining conditions for single pass 
problems which were not necessarily optimal. Hati and Rao [1976] have foimd OMC using 
this approach but have considered reduced diameter after each pass. SUMT alongwith 
DFP has been used for minimization of production cost for all passes simultaneously. The 
minimum production cost for all feasible number of passes ranging from n^m to 
obtained and the minimum out of these optimal costs is selected. Thus optima number of 
rough passes n°^ are obtained. It may be noted that this solution method can handle a 
finish pass along with n rough passes. The solution will be optimal with respect to dg only, 
and overall it is a suboptimal solution. 

The upper limit for depth of cut for the finish pass {dgmax) is obtained using the rela- 
tionship [Kals and Hijnk, 1978] 

d’snuix ~ 

where the slenderness ratio (S) is defined as ratio of width of cut to imdeformed chip 
thickness and its value is selected such that unmanageable chip flow does not occur. Starting 
from the upper limit of depth, of cut (dgmax) maximum feed (fmax) it is checked whether 
the given feed satisfies all the constraints at a given depth of cut. If all the constraints 
are not satisfied at this feed value the feed is reduced by a small percentage till all the 
constraints are satisfied. In case the constraints are not satisfied even at the minimiun 
value of feed {fmin) then the depth of cut is reduced by a small percentage and the search 
starts from the maximum value of feed {fmax}- In other words, the highest possible depth 
of cut and the highest possible feed at which all constraints can be satisfied are selected for 
the finish pass. 

Using the constant diameter model. Shin and Joo [1992] have suggested that dg should 
be kept at the least possible value dgmin- This value of dgmm gives the smallest finish pass 


1/1 



production cost. However, it does not guarantee that the overall optimal production cost 
would result through the least depth of cut selection in the finish pass. They have intro- 
duced a constant representing the ratio of tool life for rough pass to finish pass. As the 
exact value of C for multipass turning is not known before hand, introduction of such a 
constant makes the optimal solution somewhat arbitrary. They have solved an example 
taking C = 1- The safient feature of the model is that it uses tool life as the only parameter 
on which the total production cost depends. The tool life corresponding to the minimum 
production cost has been found using a single dimensional Fibonacci’s search method. The 
minimum production cost function is evaluated as sum of one finish pass and n rough passes 
production costs. The production cost for a particular pass is evaluated using the maxi- 
mum allowable feed value at that tool life by all constraints for that pass. The speed has 
been calculated using the optimum tool life and feed values. The model has been used 
to find the minimum production cost for all feasible number of passes in the range of n^m 
to Umax- Thus the optimum nixmber of passes is found. The limitations of the model are: 

(i) Selection of dsmin restricts the possibility of reducing the number of rough passes 
by selecting higher depth of cut in the finish pass. Therefore this method sometimes give 
significantly higher minimum production cost. 

(ii) The condition C = 1 restricts the selection of the optimum speed, feed, depth of cut and 
number of passes. 

(iii) The model has been made on the basis of initial diameter and therefore does not give 
accurate minimum production cost. 

(iv) All the rough passes axe of equal depth. 

Van Houtan [1981] uses a graphical search method in speed-feed space for finding single 
pass OMC. The depth of cut for rough pass has been selected on the basis of equal depth 
of cut. A CAPP system for rotational parts has been developed using the graphical search 
method [Van Houtan, 1986] but the model uses constant workpiece diameter. Hayes et 
al. [1979] have suggested a discrete variable approach to determine OMC for minimum 
production cost for n rough passes. It is based bn exhaustive search in T — n — space, 
where rit represents the number of tool changes per part. Although this method has been 
reported to be very efficient for an example with workpiece length of one meter, it may not 
be so for small workpiece lengths because the set of feasible Ut contains a large number of 
elements. 

(b) Unequal Depths of Cut Decided in Advance 

Ermer and Kromodihardjo [1981] have presented a production cost model for multipass 
turning based on constant initial workpiece diameter. Depending upon the e3q>erience, the 


15 



number of passes and the depth of cut distribution among these passes are used as input to 
the model. A GP-LP technique has been used for the determination of optimum production 
cost, speed and feed for all passes. The technique can, in principle, handle nainimization 
of a posynomial function imder posynomial constraints. Therefore, all the passes can be 
handled simultaneously. 

With the introduction of a few variables and a constraint, the problem of maximization 
of profit rate function has also been converted to minimization of a posynomial function 
subjected to posynomial constraints [Gupta et al., 1994]. The only drawback of this kind 
of OMC determination is in deciding the optimal input number of passes and depth of cut 
distribution. 

Hinduja et al. [1985] have suggested a few rules for the selection or modification of non- 
uniform depth of cut distribution in rough passes. A term "depth of cut corresponding to 
minimum production cost" has been used to guide this selection. Hinduja etal., however, did 
not use this approach in their subsequent papers and adopted equal depth of cut strategy 
[Chen et al., 1989; Arseculartene et al., 1992]. The saUent feature of the above research 
work is that the satisfactory chip-breaking region is divided into a grid and graphical search 
is carried out to deterrnine OMC. This ensures easy swarf disposal. On the other side, there 
is always a possibility of a finer grid and better solution. Further, solution obtained is not 
optimal. 

Yellowley and Gunn [1989] have presented some proof in support of two unequal rough 
passes approach. According to them two rough passes of equal depth will never give min- 
imum production cost or production time and any greater depth of cut will always lead 
to better solution. There proof is limited to two pass case only but the same should also 
hold good for any number of passes. The tool life eqn. (1.6) has been used to reach the 
conclusion and therefore the same conclusion should hold good for extended Taylor’s tool 
life eqn. (1.2) also. 

(c) Variable Depths of Cut and Dynamic Programming Approach 

Iwata et al. [1977] have used dynamic programming approach to find the optimal number 
of passes and depth of cut for each pass. They analyzed the multipass turning problem as 
a multistage decision making process. At each stage the optimal solution of the previous 
stage was used as input and a single pass solution method gave the optimum solution for 
the current stage. The advantage of this approach is that optimum number of passes are 
uniquely determined as stages required to machine the total depth of cut. The reduced 
workpiece diameter is made the input to a stage and thus varying workpiece diameter 
models can be easily, solved using DP approach. In this method the total depth of cut is 
divided into a number of small segments and the optimal distribution of depth of cut is taken 
as number of segments in each pass. The depth of cut accuracy for a pass will, therefore. 



be of the order of the dimension of the segment. Any effort to increase the accmacy would 
require large number of single pass turning solutions. For example Iwata et al. [1977] have 
used 10 divisions for a total depth of cut of 10 mm. Hence, the pass can be as accurate as 
1 mm. If the depth of cut range for a pass is 1-3 mm, it requires approximately 10 x 3 
= 30 single pass solutions. If the desired accuracy is 0.1 mm then the number of divisions 
required will be 100 and the number of single pass solutions will increase to 100 x 21 = 
2100. Any nonlinear programming method coupled with dynamic programming approach 
to solve a multipass turning problem would thus reqxiire prohibitive time. Iwata et al. have 
used SUMT along with Newton Raphson method to solve single pass problems for small 
nirmber of solutions. 

Lambert and Walvekar [1978] have developed a DP approach to solve two pass rough 
turning problems. GP was used to find single pass solutions. They found that the optimmn 
solution often contains two unequal rough passes. Unklesbey and Crieghton [1978] have 
simplified this approach and have obtained one finish pass and one rough pass solutions 
on the basis of constant workpiece diameter. The same rough pass solution is used for any 
number of rough passes. GP has been used to obtain single pass solutions. Thus the number 
of single pass solutions has been considerably reduced. The same DP approach has also 
been used for the determination of optimal number of passes for weighted sum of production 
cost and production time criteria [Agapiou, 1992 b]. A NLP method (Neldermead simplex 
method) is used to find single pass solution. This work is very similar to that of Iwata 
et al. [1977]. The same approach has also been used for obtaining the optimal machining 
conditions for face milhng operation recently [Sonmez et al., 1999]. 

(d) Variable Number of Passes and Depths of Cut 

Since the optimum total number of passes are not known in advance, hence it may be 
treated as a decision variable. Further, the depths of cut for each pass are also unknown at 
the beginrdng and therefore these are also decision variables. In this case, the total depth 
of cut constraint is added to the constraint set. 

Chua et al. [1993] have constructed a multipass turning model for rough passes on the 
basis of initial workpiece diameter for each pass. The objective function (production time 
per piece) was a polynomial function subjected to one equality constraint in the form of total 
depth of cut and two monomial constraints for each pass. A NLP method, the sequential 
quadratic programming (SQP), has been used for optimization. The problem has been 
solved in two stages. In the first stage n is treated as a continuous variable and a real 
optimum number of rough passes (n'^’’) is found on the basis of constant initial workpiece 
diameter, equal depth of cut for all passes and the same cutting conditions, i.e., the same 
speed and feed in all passes. In the second stage the same objective function is minimized for 
[j.jc^r]+ j-j.joptrj- model has been made on the basis of changed workpiece 
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diameter in each passes and with no restrictions on speed, feed and depth of cut for any 
pass. The minimum out of these two solutions is selected as the solution. This is obviously 
not an optimal solution to the problem because there is no guideline to the selection of 
finish depth of cut. 

Kee [1996] has suggested that the solution obtained by the above method has high 
chances of converging to the local minimum at the first stage when n is treated as a con- 
tinuous variable. Further, there is no evident termination point to the optimal number of 
passes considered in the search. He has suggested an alternative approach to the solution 
of multipass rough turning problem. It is based on the study of theoretical economic trends 
for number of passes and depth of cut distributions, fully utilising the constrained single 
pass optimization strategy to establish speed and feed for each pass to reach the global 
optimum The main features of this work are: (i) it contains a logical method for nar- 
rowing down the range of number of passes, i.e., rimin to to locate the global optimum 
number of rough passes (ii) it is based on calculating the production time for multipass 
turning model with changed workpiece diameter at each pass for a given and depth of 
cut distribution, and (iii) at the last stage, the depth of cut distribution is refined using a 
penalty method [Schuldt et al., 1977]. 

It has been claimed that the above work successfully locates the global optimum number of 
rough passes. The distribution of depth of cut is obtained using a penalty function method 
and it may be local minima. In the absence of guidelines for the selection of finish pass 
depth of cut, the solution obtained is not optimal. 

Tan and Creese [1995] have developed a complete model for the calculation of optimal 
total number of passes (n^) alongwith the speed, feed and depth of cut for each pass. They 
have made the production cost model on the basis of initial workpiece diameter at each pass 
for maximum possible number of passes Up^ax = (^rmox+i) and introduced the same number 
of binary [0,1] variables. A few constraints have also been added to the original problem 
with practical constraints. Although the practical constraints are in monomial form, they 
can assume any other form. The additional constraints ensure that the optimal solution 
satisfy the binary nature of newly introduced variables while these are treated continuous 
in the range of zero to one during the solution method. Sequential Linear Programming 
has been used to obtain the solution for NLP model. A feasible solution is reqxiired as the 
starting point and radius of bound or move limits for each variable are to be defined to start 
the procedtue. The selection of move limits is a trial and error process. The solution method 
although quite simple (conceptually and numerically) it is not recommended as a technique 
for general design optimization environment due to lack of robustness and uncertainty in the 
specification of move limits [Arora, 1989]. It appears to be the only complete general model 
for multipass turning based on starting diameter of the workpiece at each pass. The value of 
is obtained by counting non-zero binary variables in the optimal solution. Further, the 
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speed, feed and depth of cut associated with each pass corresponding to non-zero variable is 
also available with the solution. The model can be solved using any of NLP search methods 
viz. SUMT with DFP, SQP, etc. This method again gives local minima and has no evident 
termination point for depth of cut as in the model proposed by Chua et al. [1993]. 

(e) Largest Cut and Heuristic for Turning of Stepped Shafts 

It is an important decision to club machinable volumes horizontally, vertically, or in circular 
manner and then decide about multiple passes, during rough turning of stepped shafts. The 
production cost and production time is affected by this grouping of machinable voliunes. 
The recent available research works use a common practice of having the largest possible cut. 
[Jha, 1996; Prasad et al., 1997]. In a research work related to milling a heuristic based on 
minimmn production cost per unit volume has been used to obtain the Tninimmii production 
cost [Yellowly and Fisher, 1994]. The solution to the milling problem suggest that sometimes 
vertical or circular grouping of machinable volumes may lower the production cost. 

1.3 Scope and Organization of Thesis 

The review of existing literature reveals the following: 

1., The problem of determining the optimal machining conditions for turning as well 
as other machining operations has been drawing attention of various researchers since the 
beginning of the century. However, much of the literature deals with generalization of mod- 
els, and adoption and development of newer solution methodology. There appears to be an 
increasing interest in the determination of optimal machining conditions although some of 
the processes are quite old. '' 

2. The minimum production cost and minimimn production time criteria are most com- 
monly used for the determination of optimal machining condition. Both production cost and 
production time functions are equivalent from the point of view of optimization. The later 
can be obtained from the former by neglecting some of the costs and treating others as unity. 

3. The optimization problems are formulated on the basis of deterministic as well as 
probabilistic tool life equations and cost parameters. The deterministic problems are easier 
to solve and several deterministic model solutions claimed to yield results close to those 
obtained using stochastic models. 

4. Most of the reported works are based on two eisstimptions: (i) constant initial work- 
piece diameter and (ii) predetermined uniform/non-uniform depth of cut distribution during 



all passes. Depending on the assumptions, different sets of optimal values of speed, feed 
and depth of cut for each pass are obtained. Further, it may be cumbersome to implement 
these optimal sets of cutting conditions during manual setting of machining conditions but 
in the present environment of automated machining with NC/CNC, the setting of different 
speeds, feeds and depths of cut in various passes would not pose any difficulty. 

5. A third assumption of equal tool life in rough as well as finish pass has also been 
used by some of the researchers to simplify the solution of multipass turning problem with 
constraints. 

6. The majority of technical constraints described in section 2.3, except the constraints 
of tool wear and effective chip-breaking region, are monomial constraints containing con- 
stants and exponents for speed, feed and depth of cut. The estimates of these constants 
and exponents can be made through experiments. The determination of optimal machining 
conditions would therefore greatly depend upon the quality of these estimates. 

7. There are two t 3 q)es of approaches to find the optimal machining conditions for turn- 
ing. In the first one, the nrunber of passes, speeds, feeds and depths of cut for all passes 
are treated as decision variables and their optimal values are obtained from the solution 
of the model. In the second approach, some of the decisions regarding number of passes 
and depth of cut distribution are assumed to be known based on some practical knowledge. 
Their values are given as input to the model and remaining decision variables are obtained 
from the solution of problem. 

8. The various types of optimization methods and algorithms used to obtain optimal 
machining conditions are as follows: 

(i) In a multipass ttuning problem all the machining conditions including the number 
of passes, speed, feed and depth of cut for each peiss are treated as decision variables. The 
objective function as well as constraints are linearized about a feasible starting point and 
linear programming is sequentially applied to better the solution till the optimum is found. 
The determination of a feasible starting point requires solution of a separate optimization 
problem and application of sequential linear programming requires careful selection of move 
limits. 

(ii) The multipass turning problem with constraints is converted into unconstrained 
problem using a penalty function method. The resulting unconstrained problem is solved 
using SUMT applying any non-linear search method such as Powell’s method, Davidon 
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Fletcher Powell’s method or Neldermead Simplex method, etc. These methods require; (a) 
predetermined value of number of passes, (b) experimentation with a few parameters during 
SUMT application and (c) possibly a feasible starting point. 

(iii) In some of the works the constrained optimization methods viz. generalized reduced 
gradient (GRG) and sequential quadratic programming (SQP) have been applied for the 
solution of multipass turning problem with constraints. GRG method requires a feeisible 
starting point, whereas SQP may be started with any arbitrary point and also without any 
knowledge of parameters for optimization. Therefore, SQP method should be preferred over 
other methods. 

(iv) A special type of nonlinear programming method - geometric programming (GP) 
has been widely used for single as well as multipass turning optimization where the objective 
function and the constraints happen to be posynomial. The GP based methods may or may 
not require a feasible start point. Sometimes the knowledge of depth of cut distribution 
and number of passes is required to solve the multipass problem. 

(v) Another method for solution of multipass tinning problem is to break it into several 

single pass turning problems using predetermined number of passes and depths of cut or 
use of dynamic programming (DP) which automatically decides the optimal number of 
passes and depth of cut distribution. The resulting single pass turning problems may be 
solved using any non-linear programming method, or any other method using differential 
calculus and graphical search over effective chip breaking region defined in section 2.3. An 
assumption of equal tool life is used in the method alongwith assumptions of equal workpiece 
diameter and predetermined uniform or non-uniform depth of cut. The resulting multipass 
turning problem can be solved using single dimensional search in terms of tool fife. At 
a given tool life, the maximum feed corresponding to the constraints yields the minimum 
production cost. The method does not require a feasible point. If the three assumptions 
made above are removed, it turns out to be an easy method than those described in (i) to 
(iv). • 

1.3.1 Scope 

The mim'nmm production cost or minimum production time criteria are most commonly 
used for the determination of optimal machining conditions. Both of these criteria are 
equivalent from the point of optimization. As the later can be derived from the former 
one, the minimum production cost criterion has been adopted in the present work. The 
deterministic tool life and cost coefficients have been used to make multipass turning model. 
As mentioned earlier, there are several approaches to solve multipass turning problem. 
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These include: (i) linerization of objective function as well as constraints about a feasible 
initial point and use of linear programnaing sequentially to better the solution, (ii) conversion 
of constrained optimization problem to unconstrained optimization one and to solve the 
resulting unconstrained problem, (iii) use of constrained optimization methods starting 
from feasible or arbitrary point, (iv) transformation of complex multipass turning problem 
into several single pass turning problems using assumptions of equal workpiece diameter, 
equal tool life, predetermined number of passes and uniform or non-uniform depth of cut in 
each pass, and (v) use of DP to find the optimal number of passes and optimal distribution 
of depth of cut using single pass solutions. 

In the present work, the last two approaches have been considered as they are based 
on decomposing the complex multipass turning problem into simple single pass turning 
problem. The available single pass solution methodologies are: (i) nonlinear constrained 
or unconstrained approach, (ii) geometric programming applicable to posynomial function 
with posynomial constraints, (iii) graphical methods using differential calculus to search a 
minima over (v — d — f) space in which the objective function as well as the constraints are 
described, and (iv) determination of optimal tool life corresponding to Tm'ni'mnTri production 
cost using any one dimensional search method and evaluation of TninimiiTn production cost 
by finding the maximum permissible value of feed under the constraints. The last method 
has been used to obtain single pass solution in the present work because this method works 
efficiently under the monomial constraints. 

The number of passes may be treated as decision variables and the optimal number of 
passes may be obtained from the solution of optimization problem. This number may also 
be predetermined based on experience. Two models have been developed in the present 
work. The first one starts with the maximum nmnber of total passes and finds the optimal 
number of passes alongwith optimal subdivision of depth of cut among passes. The model 
has been developed on the assumption of equal (initial) workpiece diameter and equal tool 
life in all passes. This model yields lower minimiun production cost than that obtained 
using a commonly practised strategy in which the finish pass is made at the lowest possible 
depth of cut and all rough passes are made of equal depth. The optimal subdivision of depth 
of cut obtained using this 'model does not reveal the order in which the depths of cut for 
rough passes should be used. The second model makes use of dyna m ic programming and, 
therefore, varying workpiece diameter model can be handled easily. Further, the optimal 
number of passes are determined urdquely as stages required to machine the total depth 
of cut. This model clearly discriminates the order in which the depths of cut should be 
applied. The assumption of equal tool life in all passes has also been dispensed with. The 
solution methodology for single pass turning problem has been selected such that once a 
solution to single pass turning has been obtained for a given depth of cut, the solutions 
to problems of the same depth of cut with changed workpiece dimensions can be obtained 
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using parametric analysis. While using dynamic programming to find the optimal number 
of passes and optimal subdivision of depth of cut among passes, the single pass turning 
problems with changed workpiece diameter need not to be solved again and again as has 
been done in the past. The solution to single pass turning problem with changed workpiece 
diameter can be obtained using simple calculations. 

During turning of stepped shafts, various possible sets of machinable volumes and their 
production costs may be obtained on the basis of horizontal, vertical and circular clubbing 
of machinable volumes. The minimum production cost of a stepped shaft depends upon 
the optimal grouping of machinable volumes. There is a need to find the optimal grouping 
of these machinable volumes using a heuristic method or a mathematical model to obtain 
minimum production cost. A heuristic based on miniTniiTin production cost per uni t volvune 
has been suggested to find the optimal grouping of machinable volumes. A mathematical 
model has also been suggested to obtain optimal grouping of machinable volumes yielding 
minimum production cost. 

1.3.2 Organisation of Thesis 

In Chapter 2, the time and cost calculations have been presented to form a production 
cost ftmction for multipass turning problem. The turning operation is restricted by various 
practical constraints viz. cutting power, cutting force, etc. and limits on cutting variables. 
The functional relationships of various constraints in terms of speed, feed and depth of 
cut have also been presented in this chapter. It may be noted that majority of these 
constraints are monomials. Finally, Model 1 heis been formulated to minimize production 
cost in multipass turning with monomial practical constraints alongwith total depth of 
cut constraint for a given number of passes (n^). The problem is complex without the 
knowledge of optimal number of passes (n^) and optimal distribution of depth of cut, and 
due to varying workpiece diameter in each pass. The multipass turning problem with all 
these constraints can be solved in single stage by considering all the passes simultaneously, 
or it can be broken into various single pass turning problems using some knowledge to 
decide the number of passes and depth of cut distribution before optimization. A further 
simplification has been obtained using the assumption of equal tool life in all passes. Using 
DP, the multipass turning problem can be viewed as multistage decision making process 
where each pass constitutes an stage. DP can take into account the varying workpiece 
diameter and uniquely gives optimal number of passes and various dj for i = 1, 2 . . . . 

The main difficulty in the application of DP is that it requires several single pass solutions. 
Single pass turning or multipass turning problem may be solved alongwith constraints using 
any NLP method viz. SUMT with Newton Raphson or DFP, GP, GRG and SQP. Some 
of these methods require a feasible starting point. The determination of feasible starting 
point requires solution of a series of optimization problems. Therefore, the effort has been 
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made to select those methods which do not require feasible starting point. 

A solution methodology is presented in Chapter 3 for minimization of production cost 
for single pass turning problem with monomial constraints. The methodology is based on 
the determination of optimal tool life over working range of tool life using golden section 
method. It has been observed that the maximum allowable feed under the constraints 
for a given tool life minimizes the production cost at that tool Hfe if the depth of cut is 
fixed. In the beginning, attempt has been made to find the optimal number of passes and 
optimal subdivision of depth of cut among all passes. The assumptions of equal tool life 
and same initial workpiece diameter in all passes have been used to simplify the problem. 
A hnear integer programming model has been presented to solve the multipass turning 
problem using LINDO software. It has been observed that optimal tool life, speed, feed 
and depth of cut depend only on depth of cut under the given constraints for single pass 
turning. Fiuther, these do not depend on the workpiece diameter or length. The variable 
component of minimum production cost is directly proportional to workpiece dimensions. 
This property has been used to find the minimum cost for changed workpiece diameter or 
length and it sigmficantly reduces the required number of single pass solutions during DP 
solution of multipass tiuning problem. The second model described in the same chapter is 
based on DP. This model successfully deter min es the optimal number of passes and depth 
of cut distribution alongwith speed and feed for individual passes. 

In Chapter 4 some numerical results have been presented. The first two examples 
have been Solved using integer programming model. It has been found that the commonly 
practised strategy of using small depth of cut in the finish pass and all rough passes of 
equal depth of cut does not yield the optimum number of passes or optimmn depth of cut 
distribution. Therefore, the speeds and feeds associated with this solution are also not 
optimal, whereas integer programming model gives optimal niunber of passes and optimal 
depth of cut distribution alongwith speeds and feeds for all passes. It has been found that 
uneven distribution of depth of cut gives reduced minimum production cost. There are 
further chances to improve the solution because it is based on assumptions of equal tool life 
and same workpiece diameter in all passes. The next two examples have been solved using 
dynamic programming model. It has been found that the tool hfe at optimiun solution is 
different for finish and rough passes. This model takes into consideration the effect of varying 
diameter also and provides better results. The last example give comparison of execution 
time for single as well as multipass turning problem with other nonlinear programming 
methods viz. SUMT with DFP, GRG and SQP. 

In Chapter 5 stepped turning of rough shafts has been considered. Machinable volumes 
are formed by drawing horizontal and vertical fines firpm the ends of various cylindrical 
volumes. Clubbing of machinable volumes horizontally or vertically becomes an important 
decision to minimize the production cost. The normal practice is to cut the largest possible 
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length and therefore possibilities of clubbing volumes vertically <w- eliminated. Vertical 
clubbing of volumes may reduce the number of passes as well as minirmim production 
cost. A heuristic based on the minimum production cost per unit volume and an integer 
programming model has been suggested in order to find the optimal clubbing of machinable 
volumes. The results, in general, show that the heuristic method yields better minimum 
production cost than that obtained from common practice and the mathematical model 
gives still better solution. 

In Chapter 6 the conclusions derived from the present work have been summarized. It 
also gives scope for future research as well as possible integration of the present work in 
CAPP systems for machining. 


/ 
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Table 1.1: Some of the CAPP systems for rotational parts using OMC. 
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Table 1 .2 : A summary of literature in the field of optimization of single or multipass turning. 
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Table 1-.2 (continued) : A summary of literature in the field of optimization of single or multipass turning. 
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Table 1.2 (continued): A summary of literature in the field of optimization of single or multipass turning. 
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Table 1 .2 (continued); A summary of literature in the field of optimization of single or multipass turning. 
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Keys for Objective Functions 


MC : 

Minimum Production Cost per Piece 

MRR : 

Maximum Material Removal Rate 

MP : 

Maximum Profit Rate 

MX : 

Minimum Production Time per Piece 

MC+MT : 

Minimum Weighted Avg. of Production Cost per 
Piece and Time per Piece 

M (vf) 

Maximvun of Product of Speed and Feed 

VMR/TW 

Maximum of Material Removal Volume/Tool Wear 

GP (d+MT) 

Goal Programming (Depth of Cut and MT) 

Keys for Number of Passes 

S 

Single Pass Turning 

D 

Double Pass Turning 

M 

Multipass Turning 

Key for Model Type 


D : 

Deterministic 

P : 

Probabilistic 

Keys for Constraints 

CB 

Chip - Breaking Region 

CF 

Cutting Force 

CT 

Cutting Time 

DA 

Diametral Accuracy 

DS 

Dynamic Stability 

M 

Torque 

NPC 

Number of Pieces 

P 

Power 

SCR 

Stable Cutting Region 

SF 

Surface Finish 

T 

Tool Life 

Th 

Tool Thrust 

TH 

Temperature at Tool Chip Interface 

TW 

Tool Wear 

U 

Unconstrained Optimization 

W 

Workholding Constraint 

Keys for Solution Technique 

CGM 

Conjugate Gradient Methods 

DC 

Differential Calculus 

DC Comp 

DC and Computational 

DC Graph 

DC and Graphical Search 

DC IT 

: DC and Iterative 
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DCd-f 
DC Simu 
DC Trans 
DP-GP 
DP-GP(DC) 
DP-SUMT-NR 
DP-SUMT-NMS 
EE 

ES-SR 

LFM 

GLP 

Exhaust Search 
GP 

GP-Analytical 
GP-LFM 
GP-LP (SP) 


GP-Stochastic 

IT-Comp 

LIP-(Opt-T) 

Opt-T 

SA 

Powell 

Simu-GINO 

SEP 

SLQ 

SUMT-DFP 

SUMT-DC 

Key for Methodology 


DC and Search over d-f Diagram 
DC and Simulation 
DC and Transformation 
Dynamic Programming with GP 
Dynamic Programming with GP (using DC) 

DP using SUMT with Newton Raphson 
DP using SUMT with Neldermead Simplex 
Empirical Equation 

Expert System and Storage and Retrieval 

Lagrange - Function Method 

Goal Programming 

Exhaustive Search 

Geometric Programming 

GP with Analytical Method 

GP with Lagrange’s Function Method 

GP with Linear Programming using Separable 

Programming 

GP with stochastic Programming 

Iterative and Computational 

Linear Integer Programming using Opt-T 

Optimal Tool Life Search 

Simulated Anealing 

Powell’s Method 

Simulation using GINO Package 

Sequential Linear Programming 

Sequential Quadratic Programming 

Sequential Unconstrained Minimization Technique 

using DFP 

SUMT for Refinement of Results Obtained Search 
using DC - Search 


1 

2 

3 

4 

5 


, Single Pass Turning 

Small Depth in Finish Pass and Equal Depth in All 
Rough Passes 

Depths of Cut are Predecided for All Passes 
Depths of Cut are Determined on the Basis of Equal 
Workpiece Diameter 

Depths of Cut are Determined on the Basis of Varying 
Workpiece Diameter 


Key for Inputs Parameters 

d : Depth of Cut 

f : Feed per Rev. 

np : Nvunber of Passes 
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rvcy ui \^uLpuL voiiciuicii 


d 

Depth of Cut 

f 

Feed per Rev. 

n 

Number of Rough Passes 

np 

Total Number of Passes 

nt 

Number of Tool Changes per Part 

r 

Nose Radius 

V 

Key for Tool Life 

Speed 


Modified Golding’s Tool Life Equation 
Extended Taylor’s Tool Life Equation 
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Chapter 2 

Modelling of Turning Operations 


2.1 Introduction 

The miniinura production, cost or nnininuiTn production time criteria \s most commonly 
used for the determination of optimal machining conditions. The various constraints for a 
pass are defined in terms of decision variables viz. speed, feed and depth of cut corresponding 
to each pass. In multipass turning the only constraint connecting all the passes is total 
depth of cut. The mathematical treatment of optimization problem requires a suitable 
model and solution methodology to determine the optimal machining conditions viz. the 
optimal number of passes, speed, feed and depth of cut for each pass. For simpficity, straight 
turning operation of a cylindrical workpiece has been considered. The problem is to find 
the expression for production time and production cost and the constraints for single as well 
as multipass turning. It is assumed that, in general, one final finish pass and a few rough 
passes are required to complete a job. The corresponding constraints have to be identified 
and described. The next step will be to formulate the optimization problem and to develop 
a solution methodology. 

Section 2.2 of this chapter deals with the basic model of turning. The expressions for 
objective functions viz. production time and production cost are derived on the basis of 
constant as well as varying workpiece diameter. In section 2.3 the description about various 
constraints is given. The statements of various single as well as multipass turning problems 
have been presented in section 2.4. 

2.2 Objective Function 

This section deals with the development of production time and production cost functions 
for multipass turning of cylindrical bar of initial diameter Do and length L in rtp nximber 
of passes with last pass as the finish pass. Subscript i is used to specify the pass number 
which varies from 1 to Up. 
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2.2.1 Production Time per Piece (t^) 

The production time per piece consists of preparation time (tpi), idling time (tai), machining 
time (tmi) and tool changing time {td) for all passes. 

(i) Preparation time {tpi) is the time required to prepare for loading, unloading, etc. for 
pass. 


(ii) Idling time {tai) consists of tool approach and depart time for pass and is given 
as 

^ai — huL + /l2i (2-1) 

where hn and h 2 i are constants related to tool approach and depart times. The first term 
denotes the time required for tool travel from a set point and depends on workpiece length. 
The second term h 2 i is the time for engagement and disengagement of tool with workpiece. 
Idling time consisting of hu and h 2 i, in general, is the same for all passes if these passes are 
machined on the same machine tool. Idling time for various passes will not be the same if 
these passes are taken on diflferent machine tools. 


(iii) Machining time (tmi) is the time needed for actual cutting in pass. If the diameter 
of workpiece at the start of pass is (A-i), then tmi is expressed as 


•wDi-iL 

lOOOvifi 


( 2 . 2 ) 


where Vi and fi are speed and feed respectively for pass. The intermediate diameters 
Di-i are calculated as ^ • 

A-i = T>o - 2 £ 4 (2.3) 

l2=i 

where dj is the depth of cut removed in pass. 


(iv) Tool changing time per tool or per tool tip {td) is the time required for tool or tool 
tip replacement for pass. The tool changing time per piece {td) for pass is defined as 

= U— (2-4) 

where Ti is the tool life for cutting of pass. It is related to the speed, feed and depth of 
cut by tool life eqn. (1.6) which assumes the following form for a pass i: 

= Ctu (2.5) 

Equation (1.6) reduces to extended Taylor’s tool life eqn. (1.2) as described earlier. For a 
specific pass i it is represented as 

= Cn ( 2 . 6 ) 
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Hence the total production time tto is 

i=np i=7ip 

= E (V + + E + iei) (2.7) 

i=l i=l 

The first two terms of the above expression are constants for a given workpiece length and 
number of passes. Therefore in actual minimization of production time these may be kept 
out of objective function and may be added to it after optimization. 

2.2.2 Production Cost per Piece (Cto) 

The production cost per piece (Cto) consists of material cost, preparation cost, idling cost, 
machining cost, tool changing cost and tool cost. 

(i) Material cost (itmc) is the cost of raw material per piece. 

(ii) Preparation cost (t^i) is the cost of loading and imloading per piece for pass. 

(iii) Idling cost (Cai) is the cost of tool approach and retract per piece for 

(iv) Machining cost (Cmi) is the cost of machining per piece for pass. 

(v) Tool changing cost {C^) is the cost of tool or tool tip changing per piece for pass. 

(vi) Tool cost per piece (Cu) is the cost of tool per piece for pass. 

In addition, the overhead cost (ka) includes the direct labour and indirect overhead 
cost per minute for pass. The indirect cost is necessary to produce a workpiece and it 
comprises of depreciation cost of machine tool, general administration expenses, etc. The 
overhead cost during machining (kmi) includes the cost of cutting oil, electricity charges and 
any other charges for a specific machine. The cost of tool per edge is represented by ka- 
The various cost components (ii) to (vi) can be evaluated as 


Cpi 


(2.8) 

Cai 

= kii[hiiL + /i2i) 

(2.9) 

^mi 

= (Mi “b ^mi)^mi 

(2.10) 

Oai 


(2.11) 

Cti 


(2.12) 


The total preparation cost per piece (Cp), the idling cost per piece (Ca), the machining cost 
per piece (Cm), the tool changing cost per piece (Cc) and the tool cost per piece (Ct) for 
multipass turning are expressed as the sinn of preparation cost, idling cost, machining cost, 
tool changing cost and tool cost for individual pass, i.e., 

i=np 

C, = ECW (2.13) 

i=l 

i=np 

C. = EOs (2-14) 

i=l 
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(2.15) 


a 


i=l 

i=7ip 

EOd 

i~l 


i=Tip 

Ct = E 

i=l 

and the total production cost per piece (Cto) is 


Cto — ^mc + Cp + Ca + Cm + Cc + Cj 


(2.16) 

(2.17) 


(2.18) 


Cto ^mc + E ^“0 "I" ^^[{kli + kmi)tmi {klita + kti){tmi/Ti)] (2-19) 

i=l j=l 

The first three terms of the cost eqn. (2.19) are constants for a given workpiece length, 
number of passes and cost elements. Therefore these may be kept out during actual mini- 
mization of the production cost function and may be added later. 


2.3 Constraints 


In the real life, tiuning operation is restricted by various constraints listed below. 


2.3.1 Total Depth of Cut 


Let dti be the remaining depth of cut before a pass i. The total depth of cut to be machined 
in Up passes is represented by dti and it is equal to the sum of depths of cut removed in 
individual pass, that is. 


Tip 

d(i = E^ 


( 2 . 20 ) 


i=l 

where di represents the depth of cut removed in pass. 


2.3.2 Cutting Force 

Total cutting force (Fc) is approximated by the relationship [Chua et al., 1993] 

Fc = (2.21) 


where K[, cki, (3i and 71 are constants. The value of speeds, feeds and depths of cut for 
machining should be selected such that the value of cutting force is less than its maximxun 
value (Fcmax)- Therefore the constraint for cutting force assumes the form 


where K, = K'jF^a.- 


< 1 


( 2 . 22 ) 


38 



2.3.3 Cutting Power 

T1i 6 cutting power (i^) is related to machining conditions as [Amerago and Russel, 1966] 

P = K'^v<^^f^d^ (2.23) 

where oc 2 , p 2 y 72 are constants. If the maximum available power (Pmai) is fixed, it is 
possible to work at speeds, feeds and depths of cut which satisfy the relationship 


where K 2 = K' 2 / Pmax 


K2V‘^^f^<P^ < 1 


(2.24) 


2.3.4 Torque 

The spindle torque (M) is related to decision variables as [Hayes et al, 1979] 

M = (2.25) 

where iifg, 0:3, ^3 and 73 are constants. The constraint for the torque is expressed as 

< 1 (2.26) 

where = K'^Do/Mmax Mmax being the maximum allowable torque. 

The above inequality gives the lower bound on speed, feed and depth of cut because the 
initial diameter Dq has been used in the calculation of Kz- 


2.3.5 Tool Life 

The extended tool life equation is represented as 


T = 


Ko 


v°^f^d 


70 


(2.27) 


where 


Ko = Oq = (l/rii); A = (^aM); 70 = (^3/^1) 

It may be desirable that tool life does not fall below a certain minimum value (Tmin)-, be-, 

T > Tmin 


or 

• < 1 (2.28) 


where = Tmin/Ko . 

In case the tool is engaged for more than one operation at different cutting conditions, the 
conditions giving the least tool life may be used for restricting minimum tool life. Another 
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method of restricting tool life is to restrict the minimum number of cuts of equal depth 
(Ncmm) which can be machined in one tool life, i.e., 

Tvf 

■nDoL - ( 2 . 29 ) 


or 

'^^O^^anin ^ .. 

vfT - 

where N^min is an integer only which can assume the lowest value as one. 
diameter has been used so that the constraint will be satisfied for all passes. 


(2.30) 
The initial 


2.3.6 Tool Thrust 

The thrust on tool {Fth) is given as [Hayes et al., 1979] 




(2.31) 


where and 75 are constants. 

In order to avoid tool breakage the thrust on tool should always be less than its maximum 
allowable value {Fthmax), i-©-, 

< 1 (2.32) 


where K, = [K'jFtHma.). 


2.3.7 Tensile Stress on Rake Face 

The breaking of edge of a tool due to brittle fracture is called tool fracture. Fracture occurs 
when normal stress exceeds the ultimate stress value in uniaxial tension. The point of 
maximum tensile stress lies on the rake surface at a distance of approximately double the 
contact length of chip (Fig. 2.1). The tensile stress {at) on tool rake face is given by the 
relationship [Jang et al., 1992] 

at = (2.33) 

where and /3e are constants. 

Therefore the constraint for tool fracture due to stress on rake face is given as 

< 1 (2-34) 

where Ke = K'^/atmax'-, crtmax being the maximum allowable tensile stress on rake face. 
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Figure 2.1: Stress distribution in tool [Jang et al., 1992]. 

2.3.8 Tool Wear 

A tool is subjected to two kinds of wear (i) adhesion wear and (ii) diffusion wear. The 
wear rate of a carbide tool cutting steel due to adhesion (Lwi) is given by the relationship 
[Kannatey-Ashibu, 1985] 

= A:;i {Kn + -0^7 1 (2.35) 

where K'-j-y, 0:71, z/71 and 671 are constants and the wear rate due to diffusion is 

given by 

A^2 = [exp{-K;j{Kn''' + Kn'^' (2.36) 

where ATyj, AT^^, AT^.'^, 1^72 are constants and and the permitted flank wear-land value {L^,) is 
generally taken equal to 0.3 mm [Boothroyd, 1971]. 

The diffusion wear rate should be lower than the adhesion wear rate to prevent rapid 
deterioration of the cutting tool, i.e.,CJar>g , ei- al.V99i3 


Lw 2 — Lyji 


(2.37) 



2.3.9 Primary Shear Zone Temperature 

The primary shear zone temperature (0) should be less than its maximum allowable value 
ifimax) order to avoid plastic deformation of tool. Jang [1992] have presented the following 
relationship for estimating primary shear zone temperatxire: 


9 = 

where Ks, ag, /Sg and 7 g are constants. 

The constraint for maximum primary zone temperature is given as 

Ksv°^f^(P^ < 1 


where Kg 


= Ks/9 


max’ 


(2.38) 


(2.39) 


2.3.10 Temperature at Chip-tool Interface 

The chip-tool interface temperature (0^) dominates at higher cutting speeds. It is considered 
as another limit on cutting speed. In carbide tools the predominant mechanism of wear is 
due to crater. In a wide range of working conditions between light and severe craters, the 
expression for the chip-tool interface temperature is given as [Jang, 1992] 

03 = (2.40) 


where ag, /^g, 7 g, and K'q are constants. 

The chip-tool interface temperature should not exceed the softening temperature of tool 
{9 soften) otherwise rapid wear initiates, i.e., 

< 1 (2.41) 


where lifg — (.l^g/0s&/<eri)* 


2.3.11 Stable Cutting Region 

In general, it is desirable to have continuous chips with no built-up edge. So that steady- 
state cutting conditions prevail with no force fluctuations and associated difiiculties. The 
constraint used for this purpose is [Nigm et al., 1976] 

(2.42) 

or 

< 1 ( 2 - 43 ) 

where Kio, ckio and /3io are constants. 
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2.3.12 Chip-breaking Region 

The chip breaker performs satisfactorily in a given {d - f) region. The information about 
this region is supplied by the tool manufacturer for a tool insert with a chip breaker groove. 
A typical {d- f) region for chip breaking region and its approximation is shown in Fig 2.2. 
The minimum and maximum feed limits (/wim) and corresponding to a given 



Figure 2.2: Chip-breaking region and its approximation. 


depth of cut are obtained from the relationships 


and 


fminUm fmin 


fmaxlim frr, 


tan(j>i 
d dmit 

tan(j)2 


( 2 . 44 ) 


( 2 . 45 ) 


Therefore the acceptable feed should lie between feed limits fminUm fmaxUmj he., 

fminlim ^ ^ fmaxlim 



2.3.13 Surface Finish 


The surface finish can be expressed in terms of CLA value (i2„), or as peak-to-valley height 
(h). The relationship between cutting conditions and surface finish (i?„) as given by Bhat- 
tacharya et al. [1970] is 

= R, ( 2 . 47 ) 


where ifl2^ 0^125^2? and 712 are constants. 

If Ramax fs the maximum allowable surface roughness value, the constraint for surface 
finish is given as 

< 1 (2.48) 


where K \2 ^\ 2 / ^o^'^riax- 

A relationship for peak-to-valley height 
by Armerago and Brown [1969] is 


The peak-to-valley height of the finished 
value (hmax), i-e, 


in terms of feed / and nose radius r as presented 


f 

8 r 


(2.49) 


surface should be less than its maximum permitted 


8r 


^ hmax 


(2.50) 


or 


/.< {SThmax)'^ 


(2.51) 


2.3.14 Surface Integrity 

The smface integrity of a surface is quantified by (i) the depth of zone with compressive 
stress, (ii) the depth of plastically deformed zone, (iii) the value of residual stress and (iv) 
the compressive stress on the machined surface. Jang [1992] has given the relationship for 
the estimating the values for (i) to (iv) indicated above. 

The depth of zone with compressive stress should be less than its maximum value 

(.^czmax')i i-G., 

< I (2.52) 

where / dczmax) and ^13, As, 7i3 and 613 are constants. 

The depth of plastically deformed zone should be less than its maximum allowable value 

(A zmax) j i.C,, 

< 1 (2.53) 
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where Kn — / dp^max) and ai4, ^j4, 'Y14 and 614 are constants. 

The residual tensile stress should be less than its maximum allowable value atrmox, i.e., 

< 1 (2.54) 

where Ki^ = 

rmax ) and ifjg, ais, ^15, 7^5 and ^15 are constants. 

The residual compressive stress should be less than its maximum allowable value 

C/ltrCZ^ 7 

i.e., 

<1 (2.55) 

where iiTie = / ^cmax ) and ii'ig, q:i 6, ^le, 7i6 and Si^ are constants. 


2.3.15 Diametral Accuracy 

A workpiece with length to diameter ratio of more than 6, or a disc with small length to 
diameter ratio, or a thin-walled section pose certain restrictions dming ma chining . The 
depth of cut chosen should be such that the radial force causes very little deflection of the 
workpiece. The diametral accuracy (d^) should be less than double of Tnaximnnri deflection 
i^max) [Escikoglu et al., 1985], i.e., 

da > (2.56) 


where £i^rnca is given by 


Then 


^max 


KWFrLl 


EDt 


2K[,FrLl 

EDI 


< da 


(2.57) 


(2.58) 


Here E is the modulus of elasticity for work material and is the radial force and its 
value may be taken as equal to the cutting force Fc- The value of constant K'-^j depends on 
chucking and/or holding method such that = 0.6 for workpiece held between the chuck 
and centre, = 1.4 for workpiece held between the centres and = 2.4 for workpiece 
held in chuck. 

The constraint for diametral accuracy is written as 




(2.59) 


Kir = 


EDt 


where 


(2.60) 



For hollow workpiece {D‘^ - D'’) should be used in place of and a -, 7 = a, ^ /3i = 0^ and 
717 = 7i . 

It may be observed that all the constraints described in this section except the constraint 
of tool wear and chip-breaking region described in sections ( 2 . 3 . 8 ) and ( 2 . 3 . 12 ) respectively 
are monomials. All the constraints that have been used in literature for the optimization purposes 
have been summarized above. Several of these constraints are interdependent, e.g., cutting power and 
torque, tool thrust^^and torque and any one of them may be included in the process of optimization. 


2.4 Formulation of Optimization Problems 

In general, the production time per piece function and production cost per piece function 
are equivalent from optimization point of view. The methodology which is applicable to 
minimization of production cost may be applied to minimization of production time also. 
This methodology, however, may not be applicable to maximization of profit rate. The 
problem of minimization of production cost is presented in this section. The minimization 
of production time function can also be stated in a similar fashion. The various constraints 
applicable to the minimization of production cost or time function have already been de- 
scribed in section 2.3. The decision variables for multipass turning problem, as stated 
before, are the discrete number of passes {rip), continuous speed (u,), feed (/,) and depth of 

cut (di) for each pass i; i = 1 rip. The general formulation of multipass turning problem 

is based on varying workpiece diameter without any restrictions on depth of cut distribution 
and tool life in any pauss. 


2.4.1 Varying Workpiece Diameter in All Passes 

The workpiece diameter is reduced by an amount equal to twice the depth of cut after each 
pass. Therefore accurate time and cost results are obtained. In Model 1, the problem has 
been formulated on the basis of varying workpiece diameter but without any restrictions on 
depth of cut and tool life in any pass. 


Model 1 : Varying Workpiece Diameter without any Restrictions on Depth of 
Cut and Tool Life in any Pass 

The minimization of the production cost (eqn. 2.19) for multii)ass turning problem with rip 
passes is stated as: 
minimize 


-—r 

Cto — kmc + ^2 + tai) + X/ "h + ^'ti) (^mi/l<)] 

1=1 


i = n , 

E 

i^i 


where tmi, T; and Di-i as given by eqns. (2.2), (1.2 ) and (2.3) respectively are 

ivDi^iL 


movifi 
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A-1 


subject to the following constraints: 


K^rp 


^ ^ 

1 

Do -AE 

14^1 


(i) The total depth of cut constraint (eqn. 2.20), 

i=np 

Y.di = dti 

i=l 

(ii) The total cutting force constraint (eqn. 2.22) for all passes, 

Kuvrftdt^ < 1 

(iii) The cutting power constraint (eqn. 2.24) for all passes, 

K^ivrfM" < 1 

(iv) The spindle torque constraint (eqn. 2.26) for all passes, 

< 1 

(v) The minimum tool life constraint (eqn. 2.28) for all passes, 

< 1 

(vi) The tool thrust constraint (eqn. 2.32) for all passes, 

< 1 

(vii) The tensile stress at rake face constraint (eqn. 2.34) for all passes, 

< 1 

(viii) The tool wear constraint (eqn. 2.37) for all passes, 

Ljjj2i ^ LwXi 

(ix) The primary shear zone temperature constraint (eqn. 2.39) for all passes, 

< 1 

(x) The chip-tool interface temperature constraint (eqn. 2.41) for all passes, 

< 1 


(2.61) 

(2.62) 

(2.63) 

(2.64) 

(2.65) 

( 2 . 66 ) 

(2.67) 

( 2 . 68 ) 

(2.69) 

(2.70) 
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(xi) The stable cutting region constraint (eqn. 2.43) for all passes, 

< 1 

(xii) The surface finish constraint (eqn. 2.48) or (eqn. 2.51) for all passes, 

or 

fi < (Srhi maxY^ 

(xiii) The surface integrity constraints: 

(a) The depth of zone with compressive stress (eqn. 2.52) for all passes, 

< 1 

(b) The depth of plastically deformed zone (eqn. 2.53) for all passes, 

< 1 

(c) The residual tensile stress (eqn. 2.54) for all passes, 

< 1 

(d) The residual compressive stress (eqn. 2.55) for all passes, 

< 1 

(xiv) The diametral accuracy constraint (eqn. 2.59) for all passes, 

< 1 

(xv) The efficient chip-breaking region constraint for all passes, 


(2.71) 

(2.72) 

(2.73) 

(2.74) 

(2.75) 

(2.76) 

(2.77) 

(2.78) 


fmhilim i ^ fi — fmaxlim i 

(xvi) The. bounds on decision variables for all passes, 

mill ™ max 

fi min — fi — fi max 

dx mixi — dj dj max 


(2.79) 


(2.80) 

(2.81) 

(2.82) 


The various constraints related to all passes described by eqns. (2.61) to (2.78) can be 
represented, in general, as 

T < 1 (2-83) 
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where j = 1,2, . . . The term n, represents the total number of monomial constraints 
acting in an individual pass i. 

The multipass turning problem can be simplified by making suitable assumptions of equal 
workpiece diameter in all passes, deciding beforehand the depth of cut distribution among 
all passes, and keeping tool life to be the same for all passes. The resulting optimization 
models are given below. 

2.4.2 Equal Workpiece Diameter in All Passes 

The multipass turning problem is simplified to a great extent by taking initial or some 
other (mean or final) workpiece diameter value which remains the same in all passes. Initial 
workpiece diameter Do is more commonly used as the fixed diameter for all passes. Model 
1 can now be simplified as Model 2 given below. 


Model 2 : Constant Initial Workpiece Diameter in All Passes without any 
Restriction on Depth of Cut and Tool Life 

The objective function is same as in model 1, i.e., 
nodnimize 

i'-Up i=np 

Cto = kmc + + ^ai) + krrii)tmi + {kutd + kti){tmi/Ti)] 


i=l 


i=l 


where tmi is given by eqn. (2.2) after setting A.i = Do for all passes and Ti is given by 
eqn. (1.2). Thus, 

(2.84) 


and 


tmi — 


Ti 


lOOOui/i 


C ’"U 
T 


i. 22.^ Uli 


(2.85) 


subject to constraints (i) to (xvi) of Model 1. 

It is important to note that the objective function is now the sum of production costs 
for individual p<xsses minimized under the constraints for that individual pass, provided 
the optimal distribution of depth of cut is known. In other words, the total depth of cut 
constraint is the only constraint which connects the costs of various passes. 


2.4.3 Predetermined Depth of Cut Distribution Among All Passes 

The earlier solutions of multipass turning problem were based on predetermined number 
of passes and depth of cut distribution for all passes. A rule may have to be prescribed 
to select the depth of cut distribution at the outset satisfying the total depth constraint. 



Models 1 and 2, with known depth of cut distribution, have been modified for varying and 
constant workpiece diameter in the following subsections. 

Model 3(a) : Predetermined Depth of Cut Distribution with Varying Workpiece 
Diameter in All Passes 

The objective function in this case becomes 
minimize 

i=7ip 

Cto = ^ H- + (fetci + h^itrai/Ti}] 


1=1 


i=l 


and Di—\ are given by eqns. (2.2), (1.2 ) and (2.3), respectively. Now f„jj, Ti 
and Di-\ are obtained as 




T, 


movifi 

Jl. 

~X.. !l2.t liat 


h^i-1 

Di-x = Do % 

«L=l 

subject to constraints (ii) to (xvi). The total depth of cut constraint has been satisfied 
during the selection of depth of cut distribution. 


Model 3(b) : Predetermined Depth of Cut Distribution with Initial Workpiece 
Diameter in All Passes 

The objective function for this case can be written as 
minimize 

i=np i=7ip 

Gto ~ ^mc “t” ^ kli(tpi T tai) T ^ ^ [(Mi "b ^mi)tmi “b (Mi^ci "b Mi)(tmi/Ti)] 

i=l i=l 


where 


and 


Ti 


ttDqL 

lOOOuj/i 


x , 

C '"it 

T 

~£ iJ2.t liat. 


subject to constrcdnt (ii) to (xvi). The total depth of cut constraint is satisfied during the 
selection of depth of cut distribution. Models 3(a) or 3(b) can be interpreted as the sum of 
minimization of rip number of single pass turning models of the type 4(a) and 4(b) described 
below. 
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2.4.4 Single Pass Turning Problems 

The multiple turmng problem for a known number of passes and depth of cut distribution 

satisfying the total depth can be broken into a number of single pass turning problems 

This IS because depth of cut distribution clearly deEnes the starting workpiece diameters 

A-i corresponding to any pass i. The model given below uses varying workpiece diameters 
for each pass. 


Model 4(a) : Single Pass Turning Problem Based on Varying Workpiece Diam- 
eter 

The objective function for this problem becomes 
minimize 


Ctai - ku{tpi + + [(Alj. + krai)trai + {kuU + 


where 


tmi — 


Ti = 


ttA-iX 

movifi 

r<”i 

yyrp 


i- !i 2 L !la. 

A-i = Do— 2.^ 

i--± 


subject to constraint (ii) to (xvi) for pass. 


Model 4(b) : Single Pass Turning Problem Gased on Initial ^^orkpiece Diameter 

In this model the cost function for every pass is derived on the basis of constant initial 
workpiece diameter and the constraints for each pass are considered separately for the 
purpose of optimization. The objective function for this case is 
minimize 


Ctoi kn(tpi tai) + [{kii + + {kiitd + 


where 


tvDqL 


" movifi 


Ti 


Vym 


X M 


subject to constraints (ii) to (xvi) for pass. 
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2.4.5 Equal Tool Life in All Passes 

If the tool life in rough as well as finish pass is the same, the solution of multipass turning 
problem is further simplified. This assumption is used when all tools are changed after 
the same period. A one dimensional search method viz. golden section or Fibonacci’s 
method can be now applied to find a single value of tool life for all passes which optimizes 
production cost for a predeternained depth of cut distribution satisfying the total depth of 
cut requirement. The minimum of Cto function is evaluated using maximum allowable feed 
for every pass. The maximum feed value for a pass is obtained with the help of constraints 
for that pass only. A detailed procedure for finding out the maximum allowable feed is 
presented in section 3.2. 

Shin and Joo [1992] have suggested that the least depth of cut should be used in the last 
finish pass to get minimum production cost in the finish pass and all rough passes should 
be made of equal depth of cut. They have assumed equal workpiece diameter and equal 
tool life in all passes. 

2.4.6 Optimal Number of Passes (n^*) for Model 1 

The multipass turning Models 1 to 3 are solved for a given number of passes. The input rip 
may or may not be optimal. The methods used to guess the approximate optimal value of 
number of passes are given below. 

(i) The research work in multipass turning started with the selection of a small depth 
of cut {ds) for the finish pass and dividing the remaining depth of cut equally between the 
rough passes (n). The minimum of the objective function (production cost or production 
time) corresponding to each value of n is obtained. The minimum of these values is then se- 
lected as the minimum of objective function and the corresponding optimum value of rough 
passes or total passes (n^) is obtained. It may be noted that the solution obtained through 
this method depends upon the value of depth of cut initially selected for the last finish pass. 
There is always scope for improvement of the minimum objective function value because 
an uneven dei)th of cut can result in improved objective function [Yellowley and Gunn 1989]. 

(ii) In a slight modification to the above method for selection of optimal number of rough 
passes, it has been suggested by Shin and Joo [1992] that for the finish pass the depth of cut 
should be kept at the minimum value dsmm to obtain the least production cost in finish pass. 

(iii) The above two methods are based on treating the ntimber of rough passes outside 
the optimization problem. The multipass turning problem has to be solved for every value 
of number of rough passes ranging firom {rimin) to (rimax), treating the number of rough 
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passes n as a continuous variable. The cutting parameters viz. speed, feed and depth of 
cut for all rough passes may be kept equal. Model 2 is now solved on the basis of some dg 
value. The resulting real number of rough passes thus obtained guides in the selection 
of optimal number of passes from Model 1. Model 1 is solved for and '[n°P^’’]+ and 

without any restriction on the cutting variables for a given value of dg. This results in 
uneven distribution of depth of cut for rough passes. The finish pass depth of cut dg can 
also be made a decision variable at this stage in order to improve the solution. This part 
has not been tried in any of the research work reviewed in the present work. 

(iv) The solution of the above three methods depend on initial selection of finish pass 
depth of cut. Therefore, there is need for a method which does not require initial selection 
of depth of cut for finish pass or the optimal number of passes. 

(v) In a recent resemch work [Tan and Creesse, 1995], the multipass turning model 
(Model 1) has been redefined for assumed maximum nvunber of passes. A few binary vari- 
ables corresponding to each pass have been introduced and a few constraints into multipass 
turning problem have been added. These additional constraints ensure that even though the 
binary variables are treated as continuous, in the solution they appear only as binary. The 
optimiun number of passes are obtained indirectly by counting non-zero binary variables 
corresponding to each pass. Another important difference between this problem formulation 
and Model 1 is that the order of passes is reverse, i.e., i = 1 implies the finish pass and 
i — 2 to i = n'p imply the remaining rough passes. The number rip is selected such that it 
is greater than {rimax + !)■ The solution of this problem contains the optimal total number 
of passes n^^ The sequential Unear programming has been used to find the solution and 
therefore, any NLP search technique can be used to obtain the solution. The formulation 
of the general multipass turning problem as proposed by Tan and Creesse [1995] is given 
below. 


Model 5 : Multipass Turning Problem with Varying Workpiece Diameter and 
without any Restrictions — Assumed Maximum Number of Passes 


The objective function for this model can be written as 
minimize 

Cto = kmc + X/ +kmi)tmi + + Mi) 

i—1 i~l 


Ci = 7r 


Dj+iL 

lOOOvifi 


( 2 . 86 ) 


(2.87) 


where 



and 


te=i 

( 2 . 88 ) 

subject to constr^xints (i) to (xvi) for each pass i; i = i, . . . , n;. Here Df is the final diame- 
ter of the workpiece. The model contains additional binary [0,1] variables corresponding 
to each pass. First pass or finish pass is denoted by i = 1. The last rough pass is denoted 
by i = 2 and so on. The additional constraints are listed below. 


The finish pass always exists, i.e., 

”^1 = 1 (2.89) 

Further, if pass exists then (i - 1)*'‘ pass should necessarily exist, i.e., 

mi„i - mi > 0 (2.90) 

The following constraint ensures that variables rm are binary even if rrii are treated as 
continuous during optimization, i.e., 

. / 

wiidj = (2.91) 

1=1 

The zero value of binixry vjnriable corresponds to non-existence of that pass. Mathematically, 
if mi = 0 then di = 0 and if rrii — 1 then di > 0. These conditions can be combined and 
written as 

(di/100) - mi < 0 (2.92) 

where i = 1, . . . ,ny. Here the depth of cut has been divided by one arbitrary constant 100, 
assuming the maximum depth of cut in a pass will be less than 100 mm. The value of this 
constant may be chosen depending upon requirement. 

(vi) The complex multipass turning problem is broken into several single pass turning 
problems. The total depth of cut is divided into a number of equal segments. The depth of 
cut for each pass is an integer multiple of predetermined segments. The total production 
cost for multipass turning will be the sum of production cost for all passes. If the multipass 
turning problem involves Up passes, they are represented as rip stages of decision making 
process. Each single stage problem can be solved using any single pass turning method. 

The multipass or mxiltistage turning problem is shown in Figure 2.3. Each stage is 
represented by a cutting pass. The first (Up - 1) stages represent rough passes varying from 
1 to [rip — 1) and the last stage rip represents the finish pass, dn is defined as a state variable 
representing the remaining depth of cut before stage i. In other words, dti represents the 



total depth of cut to be removed and <^<( 71 ^ 4 . 1 ) is equal to zero since the total depth of cut 
been removed in rip passes. The decision variables are speed, feed and depth of cut for all 
passes alongwith the number of passes. The minimization of production cost for multipass 
turning operation can be stated as multistage decision making process as indicated below. 



^KVl) 


Figure 2.3: Multip<iss turning described using dynamic programming. 


Model 6 : Multipass Turning using Multistage Decision Making Process or 
Dynamic Programming 


For this model the objective function is 
minimize 

Cto 


i=7ip 

i=l 


(2.93) 


subject to the state transformations 

dt{i+i) = d<(i) - di (2-94) 

and constraints (ii) to (xvi) for each pass i; i = 1 , . . . , Up. 

The total depth of cut constraint has been incorporated in multistage decision making 
process by the introduction of transformation equations. The production cost Ci+\{dt{i+i)) 
at any stage (i j- 1 ) is defined as the minimmn production cost for stage (i + 1 ) to stage rip 
and it is described as 


='V 


Ci j-l 


minimum [ V) Ctm] 

l=i+l 


(2.95) 


r 







Further, 


Ciidti) = minimum + a+i(cit(i+i))] (2.96) 

subject to the constraint (ii) to (xvi). 

Now the total depth of cut is divided into a number of equal segments. The process of 
calculation starts from the last finish pass n^. The cost C,^{dt^) is calculated for various 
possible values of d(q)th for finish pass using any single pass solution methodology. The 
past practic.e has to optimize these single stage problem each time it is formed [Iwata 
et al., 1977, Agapiou, Therefore, any effort* to increase the number of segments of 

depth of cut increases the number of single pass solutions. For example, Iwata et al. [1977] 
have selected 10 divisions for a total depth of cut equal to 10 mm. As mentioned before, 
the distribution of passes can be as accurate as 1 mm if the depth of cut range is 1 mm to 
3 mm and it requires approximately 10 x 3 = 30 single pass solutions. If the depth of cut 
distribution is required to be as accurate as 0.1 mm, then 100 segments and approximately 
100 X 21 = 2100 single pass solutions will be required. Any non-linear programming ap- 
proach coupled with dynamic programming to solve multipass turning problem would thus 
require prohibitive time. 

All the models iiumtiomHi above are modified according to actual cutting environment 
such as the machine tool, cutting tool, practical constraints, bound on decision variables, 
etc. These may be the same or different for rough passes. However, in the present work it has 
been assumed that all passes are macliined on the same machine tool and all the rough passes 
are machined with the same type of tool. Although a general solution methodology should 
consider all these types of constraints, only monomial constraints have been considered 
alongwith effective chip-breaking region constraint for easy swarf disposal. In the next 
chapter two models have been presented which give the depth of cut distribution as well as 
the optimum number of passes. The first model is based on the assumption of equal diameter 
and equal tool life in all passes and has been formulated as Hnear integer programming 
problem, whereas second is based on dynamic programming. The single pass solution 
methodology has l)eon selected in such a way that single pass problems formed at different 
stages are not solved each time, instead their solution is obtained using solutions obtained 
at earlier stages. 



Chapter 3 

Solution Methodology - Multipass 
Turning 


3.1 Introduction 

The various multipass turning problems and associated models have been stated and dis- 
cussed in Chapter 2. This chapter presents various solution methodologies for multipass 
turning problems. These have been catagorized on the basis of (i) optimization method, (ii) 
treatment of number of passes as known parameter or decision variables, (iii) simultaneous 
treatment of number of passes or synthesising multipass solution using single pass solutions 
and (iv) remarks. 

Table 3.1 contains some of the popular methodologies. The first methodology represents 
commonly used strategy of using small depth in the finish pass and equal depth in all rough 
passes. The multipass problem is then solved repeatedly for various feasible values of 
number of passes to obtain the optimal number of passes. The second methodology aims 
at automatic selection of number of rough passes by treating them continuous variable 
[Chua et al, 1993]. The third methodology redefines the multipass turning problem for 
assumed maximum number of passes. A binary variable corresponding to each pass has been 
introduced and a few constraints have been added to the problem. These constraints ensure 
that binary variables may be treated as continuous during optimization and in solution they 
appear as binary only [Tan and Creesse, 1995]. All the passes are optimized simultaneously 
using sequential linear programming (SLP) method. The optimal number of passes are 
obtained by counting non-zero binary variables in the solution. It appears that there is a 
need for automatic selection of optimal number of passes as well as depth of cut distribution. 
The fourth methodology is an outcome of the present work. The production costs for a series 
of feasible depths of cut have been obtained for finish as well as all rough passes. These costs 
have been derived on the basis of constant initial workpiece diameter and equal tool life in 
all passes. The multipass turning problem has been reformulated as integer programming 
problem for assumed mayimiiTin number of passes using production costs thus obtained. A 
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few binary variables and constraints have been added to the new problem. The resulting 
problem is elaborately described in section 3.4.1 and is solved using LINDO software. The 
optimal number of passes and associated depth of cut distribution is obtained by counting 
non-zero binary variables in solution to the integer programming problem. 

In the previous chapter the general multipass turning problem has been formulated as 
multistage decision making model (section 2.4.6). The advantage of multistage decision 
making or dynamic programming approach is that it uniquely gives the optimal number of 
passes alongwith the depth of cut distribution required to remove the total depth [Iwata et 
al., 1977; Agapiou, 1992b]. This approach could not gain popularity because a large number 
of single pass solutions are required to obtain a multipass solution. This difficulty has been 
removed in the present work by obtaining solutions for single as well as rough passes for 
a series of depth of cut only once. At later stages single pass solutions are obtained using 
simple parametric analysis of the solutions obtained for the first rough pass. Thus single 
pass problem at the next stage need not be resolved. The last methodology of using mixed 
variable optimization (discrete number of passes and continuous speed, feed and depth of 
cut) has not been found in the literature reviewed for multipass turning. 

Basically there are two types of approaches to solve multipass turning problems. The first 
approach optimizes all the passes simultaneously, whereas the second approach synthesizes 
optimal multipass turning solution using single pass solutions. The present work aims at 
exploring the later approach. Section 3.2 lists various single pass solution methodologies 
and describes the selection of a suitable methodology. This is based on selection of optimal 
tool life corresponding to minimum production cost at a given depth of cut. In section 
3.3, a parametric analysis of single pass solution methodology is presented. This analysis 
helps in finding out minimum production cost for changing workpiece diameter and length 
for a given depth of cut and constraints without resolving single pass turning problem. 
The integer programming formulation of multipass turning problem is described in section 
3.4.1. Section 3.4.2 describes how the parametric analysis helps in significant reduction in 
computational efforts for finding single pass solutions during dynamic progra mmi ng based 
solution methodology of multipass problems. Three non-linear programming techniques 
viz. SUMT with DFP, GRG and SQP have been discussed in section 3.5. These shall be 
compared with DP alongwith single pass methodology in section 4.4. 

3.2 Single Pass Turning Methodology 

The various solution methodologies available for single pass turning problem are: 

(i) Conversion of constrained problem to unconstrained one using penalty function and 
sequentially minimizing unconstrained problem employing any of the unconstrained meth- 
ods. 



(ii) Solution of constrained problem itself using constrained methods viz. generaUzed 
reduced gradient method, sequential linear programming method, sequential quadratic pro- 
gramming method, etc. 

(iii) Graphical search over objective function and constraints over (v-d-f) space using 
differential calculus and some logic based on constants and exponents of tool life equation 
and constraints. 

(iv) Use of geometric progranaming based methods which started around 1970. These 
are still preferred due to posynomial objective function and monomial constraints [Narang 
and Fisher, 1993; Jha, 1996; Prasad et al., 1997; Sonmez et al., 1999]. 

(v) Classical application of differential calculus where optimal tool life is obtained and 
the corresponding speed is found for a fixed feed and depth of cut in a single pass turning. 

In constrained turning problem, where all the constrmnts are monomial, the optimal tool 
life is obtained using one dimensional search. During one dimensional search the minimnTn 
production cost is found using maximum allowable feed under the constraints at a given 
tool life and depth of cut. Once optimal tool life and corresponding maximum feed have 
been obtained, the optimal speed is found using tool life equation. 

The last method has been used to obtain minimum production cost because it does 
not require differentiation, expensive computation viz. matrix inversion or solution of non- 
linear equations, etc. Further it does not require determination of initial feasible starting 
point which may be a necessary requirement for some of the methods described in (i) to (iv). 
Therefore this method has been adopted in the present work and the results discussed in the 
next chapter will justify its selection. The single pass turning model described in section 
2.4.4 with monomial constraints and effective chip-breaking region constraint alongwith 
bounds on variables is given as: 

Cta = {tpi + iai) + KOutmi + (3.1) 

where 

_ ttDoL 
“ lOOOvifi 
J. 
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X mo 

KOi^i — kli ~t“ kmi 

K02,i — hitci + kti 

subject to: (i) ric number of monomial constraints, 



(ii) effective chip-breaking constraint 


j'minlim i — f i ^ fmaxlim i 


(iii) bounds on decision variables, 


Vi 


^ Vi < Vi 


fi min — fi ^ fi ■ 
dt min dj dj 

Here for a given pass, j = 1 . . . ric- 

Using extended tool life eqn. (1.2), u; is defined as 

Ct 


and the objective function (3.1) is redefined as 

where 


(3.2) 


C,oi = A0ud;‘‘fr^'-"^^ + A02i 

(3.3) 


(3.4) 

A02i = K0ii{tpi + tea) 

(3.6) 


It may be noted that AOji is a constant at a given tool life Ti and the exponents (1 — n 2 i) 
and nzi are always positive. Therefore the production cost for a pass is minimized by using 
the minimum permissible value of depth of cut and maximum permissible value of feed by 
the constraints, provided the tool life remains same. 

Using the above clue the minimum of objective function eqn. (3.1) is determined at 
various test points along tool life between the minimum and maximum tool life limits 
{Ti min) and (Tj ^ax)- These Umits are obtained using the relationships 
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Using aoi = 4; h 


(1.27). 
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sifand Koi = the tool life has also been defined as eqn. 
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(3.8) 


Equations (3.5) and (3.6) can now be rewritten as 
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(3.9) 

(3.10) 


Any one dimensional search method viz. golden section or Fibonacci’s method can be used 
for selection of test points. The minimum of minimum objective function value at test points 
is selected as the optimum production cost at a given depth of cut. The corresponding 
optimum cutting parameters for that depth of cut di are obtained using the steps given 
below. 


Step 1 

The maximum permissible feed for a given constraint (/t) is obtained from 

substituting for Uj from eqn. (3.7), 

(/,y = (3-12) 

KS‘Ki,<ir‘ 

where expl = "/ji “ 7 oi |^5 exp2 = /3)i ~ The exp2 is nonzero quantity otherwise the 

constraint does not limit the permissible value of feed. 

Step 2 

The maximum permissible feed for pass {fi) is obtained as 


/.* = min {f]i) for all j 


Step 3 

The minimum cost for i>ass at the test point is obtained from 

= AOiidf^ (/;)"^'“^^ + ^02i 


(3.13) 



step 4 


At a given test point Tj the cutting speed u? corresponding to minimum production cost is 
obtained by the following relationship: 


Step 5 




'■Oi 




(3.14) 


The optimum tool life for pass (T/) is located as the result of one dimensional search. 
The corresponding maximum permissible feed value and optimum cost are also 
obtained. 


3.3 Parametric Analysis 

The multipass turning model is accurately defined using varying workpiece diameter at each 
pass instead of constant initial workpiece diameter. Therefore, the solution to single pass 
turning problem with varying diameter for a given depth of cut and cutting environment is 
repeatedly rcciuired. The minimization of production cost eqn. (3.1) under the constraints 
has been obtained l)y s<decting optimal tool life using one dimensional search method. 
The minimum production cost at a test point % is obtained by selecting the maximum 
permissible feed under the constraints. The optimal cutting speed corresponding to a given 
depth of cut is found from the consideration of optimal tool life at a given depth of 
cut. Thus, minimum production cost at a given depth of cut is function of T)* only. This 
optimal tool life T/ remains the same for changed workpiece dimensions. The practice has 
so far been to formulate a new problem each time the workpiece dimensions are changed 
and solve the new problem [Iwata et al., 1977; Agapiou, 1992b; Jha, 1996]. The cost eqn. 
(3.3) is given as 




(3.15) 

where 

K02i. 

* KOii 

(3.16) 

and 

A02i — K0\i{tpi tai) 

(3.17) 

A02i is 

independent of the decision variables and tai = h^L + /i 2 i- Therefore, 



A02i = KOuitpi + huL + h2i) 

(3.18) 


During optimization this constant can be evaluated separately and only the first term of 
eqn. (3.3) is to be minimized. 
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3.3.1 Single Pass Turning with Varying Workpiece Dimensions 


The multipass turning model is accurately defined using varying workpiece diameters for 
each pass. The objective function in Model 4(a) therefore makes use of varying workpiece 
diameter instead of original workpiece diameter. It has been a practice to solve the single 
pass turning problem separately for all intermediate workpiece diameters even when the 
depth of cut remains the same in these cuts [Iwata et al., 1977; Agapiou, 1992b]. 

In this section a simple interpretation of optimum production cost obtained in the last 
section is presented. It may be noted that optimum tool life and hence optimum cutting 
parameters v- and f; do not change with changed workpiece dimensions (diameter and 
length) for the same depth of cut and imder the same cutting environment. Further, the 
new optimum production cost can be obtained using the previous results. The production 
cost given as 

C,r = + AO,, (3.19) 


may be redefined as 
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It may be noted that effort has always been to find // under constraints at a given test 
point to minimize the production cost. Constant AOI^ or A0"i only depend on the tool life 
at test point Tu Further, these are independent of the workpiece dimensions. The optimum 
production cost for subsequent pass or changed workpiece dimensions can be obtained using 
eqn. (3.20) or eqn. (3.21) and the results of previous pass, provided depth of cut and cutting 
environment remain the same. The value of constant A02i is calculated using eqn. (3.18) 
when change, in workpiece length takes place. In actual application the first term of eqn. 
(3.19) is a variable cost. The maximum allowable feed // under the constraints for a given 
depth of cut is also constant. Therefore, the optimal variable cost i corresponding to 
given depth of cut is 

i = (3.22) 

The variable cost per unit diameter corresponding to given depth of cut is 


^tovar i 
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and the variable cost per unit depth and unit length is given as 
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(3.23) 
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These conclusions are valid for any test point and are, therefore, also apphcable to the 
optimal tool life in pass i. In that case the constants AOh and A02i are to be calculated 
using r/. 

3.4 Multipass Turning Models 

The multipass turning model described in this section with monomial constraints and ef- 
fective chip-breaking region constraint alongwith bounds on variables is given as 


minimize 


2=7lp 

Cto = kffic + + ^m) + ^ KOiitmi + ^ K02i 


l—7lp 

E 

i=l 


where 


mOvifi 

T, - 


^ X !iZv 23.f 

iifOij = kii kjYii 
^2,i ~ -(- kti 


subject to: 

(i) total depth of cut constraint, 


dj = dti 


i=l 


(ii) ric number of monomial constraints, 

(iii) effective chip-breaking constraint, 

fminlim i fi fmaxlim i 

(iv) bounds on decision variables 

Vi^in ^ — ‘^i^n ax 

fimin ^ fi fi^run: 
di^in ^ di < di^nax 


1=1 


^mi 


{ 3 . 25 ) 
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where i = 1 . . - Wp and j = 1 . . . 

The multipass turning problem, in general, has the following design variables: (i) the opti- 
mum number of passes (nf ), (ii) the optimal depth of cut distribution among the passes 
and (iii) speeds and feeds for all passes. In order to find the solution to the problem, the 
assumptions viz. constant initial workpiece diameter, predetermined depth of cut distribu- 
tion, and equal tool life for all passes are used. These assumptions make the solution of the 
problem easy but at the cost of accuracy. In the present work, an attempt has been made 
to develop two models in which these assumptions have been removed in steps. An integer 
programnnng model is proposed that is based on constant initial workpiece diameter and 
equal tool life. The second model is based on dynamic programming approach without the 
above assumptions. 

3.4.1 An Integer Programming Model 

This model is based on constant initial workpiece diameter in all passes. It has also been 
assumed that all the rough as well as finish pass tools are changed simultaneously. Also, 
instead of using a rule to define the depth of cut distribution satisfying total depth of cut 
requirement, an optimiim depth of cut distribution alongwith optimal number of passes has 
been found using integer programming. 

The total production cost minimization is achieved in two phases. In the first phase, 
separate minimum costs for individual rough pass and finish pass are determined and tab- 
ulated for various fixed values of depth of cut. The values of depth of cut are selected 
from a series of depths. In the second phase, optimal subdivision of depth of cut for rough 
passes and finish pass, optimal number of passes and minimtun total production cost are 
determined using an integer programming model. 

Phase 1 

This phase consists of determining costs for individual finish or rough pass considering var- 
ious fixed values of depth of cut. A series for depth of cut is defined. The element of 
the series, for pass, di^ is ckq = dj ,„i„ + (dj max ~ d, ■TOm)(^); 9 = 0, 1, 2...mi, where rrii is 
a suitable integer. The vahie of ruj is selected considering dj min? max th.6 IcSlSt CO lint 

available on the machine tool for adjusting the depth of cut. 

Phase 2 

In this phase, an oi>timal combination of depths of cut d^‘ and d^ , i = l,2...n, optimum 
number of rough passes and total minimum production cost are determined. Ini 
tially a large value of number of rough passes is taken. The number of passes needed to 
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remove the total depth of material will be {n + 1), including the finish pass. The integer 
variables Xiq, i 0, 1 , (j 1, 2..7Tij are defined such that 
i = 0 implies finish pass 
i = 1, 2..n; implies i*'* rough pass 
^|r = 1, 2...mi; implies correspondence to depth of cut 
1 If (iiq value of depth of cut is selected in pass, 

== 0 if diq value of depth of cut is not selected in pass. 

Minimization of total production cost Cto may be stated as the following integer pro- 
gramming problem: 
minimize 

n 

Cto = Y. (3.26) 

i=0 

Alternatively, the problem can be stated as 
minimize 

n Ytii 

(3.27) 

i — 0 ^—0 

where corresponding to a particular dj,, Ciq = Cta for i = 0, l,2..n. 
subject to the constraint.s 


mo 

^ Ao, = Ifor i = l,2....n 

(3.28) 

i 0 


mi 

(3.29) 

q^O 


n mi 

(3.30) 


i=0 


Constraint eqn. (3.28) implies that there is only one do, selection for the finish pass and the 
finish pass must always be selected. Constraint eqn. (3.29) implies that there is only one 
diq selection in cjise a rough pass is selected. Constraint eqn. (3.30) implies that the sum 
of individual depths of cut are equal to total depth of cut. The above problem is an integer 
program which is solved using standard LINDO software. Two example are presented in 
the next chapter. 

3.4.2 A Dynamic Programming Model 

The previous model is based on two assumptions viz. equal workpiece diameter and equal 
tool life in all passes. The optimal number of passes are determined by coxmting nonzero 
binary variables and the optimal depth of cut corresponding to nonzero binary variable is 
determined index of the binary variable. The integer programming model is, however, 
silent on the order in which these deptlis of cut should be applied. The dynamic program- 
Diing (DP) approach is capable of dispensing with the two assumptions made in integer 
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programming model and it clearly determines the order in which the depth of 


be applied. 


cut should 


The comp ex mu ipaes tmmng problem is broken into several single pass turning prob- 
lems. The total depth of cut is divided into a number of ecual segments. The depth of 
cut for each pass is an integer multiple number of previously defined segments The total 
production cost for multipass tmning will be the sum of production costs for all passes Let 
the multipass turning problem involve at most n, passes which are represented as n stages 
of decision making process. Each single stage problem can be solved using single pals turn 
ing model described in section 3.2 and parametric analysis. The multipass or multistage 
turning problem is described in Fig. 2.3. Model 6 as stated in Chapter 2 is modified for 
monomial constraints and effective chip-breaking constraint as given below 

Each stage is represented by a cutting pass. The first (n, - 1) stages represent rough 
pass varying from 1 to (n,- 1) and the last stage n, represents the finish pass, d,, is defined 
as state variable representing the depth of cnt remaining before stage i. In other words 
do represents the total ,i,.pth of cut to be removed and is equal to zero (the total 

depth of cut has be<ui removrsl in rtp passes). The decision variables are speed, feed and 
depth of cut for all passes alougwith the number of passes). Minimization of production 


cost for nmltipass turulug oporatioa can be stated as multistage decision maVing 
the following manner: 

process in 

minmu'ze 

i=np 

Cio ^mc ^ ^ C'tcn 

i=l 

(3.31) 

subject to: 


(i) state transformations, 

1) ~ 

(3.32) 

(ii) efficient chip-breaking region constraint, 



fminlim i ^ ft ^ fmaxlim i 

(3.33) 

(iii) ric monomial constraints related to each pass. 



Kjivrfi%ji < 1 

(3.34) 

(iv) bounds on variables 

mill ^ ^ max 

(3.35) 


ft mm ^ fi ^ fi max 

(3.36) 


di mm ^ ^ max 

(3.37). 

where * = 1, . . . , n,, and j = 1, 

. . . 
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The total depth of cut constraint has been incnr^nrof ^ • 

. process by the introduction of transformation equations Th ^ decision making 

a. any stage (.• + 1 ) is defined as the nfinin^nns production 

anH it is described as ° \ r ■‘•j 


Further, 


^i+iUk(i+i)) = minimum Ctoi] 

i=i-f 1 


to stage rty 
(3.38) 


(3.39) 


C,(dn)= minimum [Coi + Ci+,(d,(,^„)] 

subject to constraint eqns. (3.33) to (3.37). 

The total depth of cut is divided into number of equal segments. The process of calcu 
lation starts from the last finish pass n Cost C tw ^ • i i ^ 

^lu. of depth for fimsh pass ^ing the single pass solution methodology described in sec- 
ton The vanous cost coefheients and are calculated. Now for the first rough 
pass the baeJeward recursive relationship eqn. (3.39) is used to find Cm - l)(i, _ 1 aL 
constants „ and are calculated for various deptL of cutfor mugh 

p^s Once these constants „ and have been obtained for first rough pass 

these remmn the same for all rough passes for the same depth of cut and cutting environ- 
ment. Constant „ can be evaluated using parametric analj^is for any single stage 

problem formed during use of recursive relationship (3.39) for stage (n, - 2) to 1. The past 

practice has been to optimi 7 ,e a single stage problem each time it is formed [Iwata et al 
1977; Agapiou, 1992b). 


3.5 Multipass Turning Solution Methodologi 


les 


The multipassss turning model described in this section with monomial constraints and 
ffective chip brejilcing region constraint alongwith bounds on variable is given as 
minimize 


where 


7V 

1=1 

TrDi-iL 

" lOOOtii/i 

-L 
nn 

~L 

KOii = kii + kfrii 

K02i ^ kijtd 4 “ kii 


f i=np 

CuA = K, + 53 ku{t^ + t„i) + 53 + 53 iTOzi 

» I i=i i=i 


Ti 


Ti 


(3.40) 


ao 



subject to; 

(i) total depth of cut constraint, 

^ di — = 0 

i:l 

(ii) ric number of monomial (;onstraints, 

< 1 

(iii) effective chip-l)reaking constraint, 


(iv) bounds on (hadsion variables 


where i = 1 . . . rip and j = I . . . 71^. 


fmmlim i ^ fi 

fi — fmaxlim i 


min ^ 


max 


fi min. ^ fi 
fi ^ fi max 
d» mill (h 
dj ^ (ii max 


In the previous st'ct.ion a dynamic programming based solution methodology for multi- 
pass turning has bemi pr<\s<mt(’id and the optimal number of passes and depth of cut have 
been obtained. The presiuxt section describes various methods other than this, used for 
obtaining sohition of nmltipass tinning problem. An attempt has been made to compare 
the computational time for the methodology described in the last section with various other 
methods presmited in this sindion. Various methods for minimizing multipass turning pro- 
duction cost arc basi'd on: 

(i) Conversion of const, rained multipass turning problem to unconstrained one using penalty 
function method and secpK'iitially minimizing the resulting problem using Davidon, Fletcher 
and Powell’s method (variable metric method). 

(ii) Solution of constrained multipass turning model using constrained methods. Two such 
methods have* bi'cn si'h^ctf'd. The first method is generalized reduced gradient (GRG) 
method which riHiuircs a ftuisible starting point. The second method is sequential quadratic 
programming (SCJP) metluxl which can be started with any arbitrary point. 



3.5.1 SotjuGntia.! UnconstrainGd IVIimmizatioii TechniquG (SUIVIT) 

The constrained multipass turning operation model described in the last section is trans- 
formed into unconstrained optimization problem using penalty function given in Appendix 
A. The resulting problem is sequentially minimized with the help of modified parameters, 
used to construct the penalty function during each iteration. Davidon, Fletcher and Pow- 
ell’s (DFP) Method has been used for unconstrained minimization. 

Davidon, Fletcher and Powell’s (DFP) Method 

This method is also called variable metric method. It can be thought as quasi Newton’s 
method and also as conjugate direction method. It is a powerful general purpose optimiza- 
tion method making use of currently available derivatives. It starts with an initial point 
(which has to be a feasible point when using interior penalty method) and a positive sym- 
metric matrix which is initially an identity matrix. The search directions are computed with 
the help of gradient of the function at the current point and optimal length is located using 
golden section method. The point thus obtained is tested for optimality. If this happens 
to be optimal the iterative search process is terminated otherwise matrix is updated and 
search is carried out using one dimensional method. The method converges quadratically 
because it is a conjugate gradient method. 

Sometimes it is not possible to find an initial feasible point easily. A separate optimiza- 
tion problem is then formulated to obtain the initial feasible solution. This initial solution 
is taken as a starting point and penalty function is formed with the help of suitable parame- 
ters. The penalty function is now minimized using DFP. The parameters of penalty function 
are changed sequentially and new penalty function is formed till no significant change in 
objective function occurs. The complete methodology is called Sequential Unconstrained 
Minimization Technique (SUMT) with DFP. This methodology has been implemented in 
FORTRAN 77 and is discussed in Appendix A. 

3.5.2 Constrained Programming MGthods 

Constrained programming methods or primal methods are those methods where a con- 
strained optimization problem is solved directly. It has been reported that one such con- 
strained method called generalized reduced gradient (GRG) requires less execution time 
than SUMT using DFP for single pass milling operation with non-linear objective function 
and monomial constraints [Eskcioglu et al., 1992]. However this method also requires an 
initial feasible point. Another constrained method which may be started from an arbitrary 
point is called sequential quadratic programming method (SQP). In this method the ob 
jective function as well as constraints are linearized about a selected point and linearized 



objective function is made quadratic by adding a quadratic step length. Thus the resulting 
quadratic programming problems are sequentially solved. 

Generalized Reduced Gradient Method 

The method hereafter abbreviated as GRG is due to Abadie and Carpentier [1969]. This 
method is a natural extension of Wolfe s reduced gradient method. Wolfe solves problems 
with linear constraints and non-linear objective function. GRG as an extension to this 
method is concerned with non-linear equality constraints. The non-linear inequality con- 
straints are converted to equality constraints by adding slack variables. The complete code 
has been prepared and implemented in FORTRAN 77. The general idea about the GRG 
algorithm and necessary details for its implementation are given in Appendix B. 


Sequential Quadratic Programming Method 

In this method the objective function as well as constraints are linearized about some point. 
The resulting linear objective function is made quadratic by addition of quadratic step size. 
The quadratic step size i.s added to eliminate the need for approximation of linear move 
limits. Thus the resulting (piadratic programming problem has been solved using augmented 
Lagrange multiplier in(d,hod. The solution of one quadratic programming problem is used to 
generate next quadratic programming problem and therefore the method is called sequential 
quadratic programming m<d,hod (SQP). More details have been provided in Appendix C. 




Table 3.1: Various methodologies for obtaining optimal machining conditions for multipass turning operations. 



Chapter 4 


Numerical Results for Multipass 
Straight Turning 

4.1 Introduction 

In section 3.2 of previous chapter a methodology to minimize the production cost for mul- 
tipass turning problem has been presented. It is based on the determination of optimal 
tool life corresponding to minimum production cost tmder monomial constraints iising one 
dimensional golden section search method. It has been observed through eqn. (3.3) that the 
minimum production cost at a given tool life for a given depth of cut is obtained by uRing 
maximum allowable feed under the constraints for single pass turning. The present work 
aims at solving multipass turning problem by using single pass solutions. This means finding 
the optimal number of passes and depth of cut distribution corresponding to all passes. The 
speeds and feeds for all passes are also to be determined using single pass solutions. Two 
models based on integer progra mmi ng and dynamic programming have been presented in 
sections 3.4.1 and 3.4.2, respectively. These models help in synthesis of results of multipass 
turning problem through single pass solutions. The integer programming model is based on 
constant initial workpiece diameter in all passes, whereas the dynamic progra mmi ng model 
is able to treat varying workpiece diameter. The minimum production cost corresponding 
to a series of feasible depths of cut for rough as well as finish passes has to be obtained for 
both the models using single pass solution methodology described in section 3.2. 

Integer programming model of multipass turning problem is formulated on the basis of 
one final finish and assumed maximum number of rough passes. Series of possible depths 
of cut for finish as well as rough passes are made. The minimum production cost for each 
depth of cut in series is obtained using the single pass solution methodology (section 3.2) and 
a binary variable is associated with each depth. The multipass cost minimization model is 
transformed to linear integer prograxnming model where the cost function is obtained using 
Diinimum production costs and variables and constraints are constructed such that only 
one finish pass and a number of rough passes are selected in order to cut the total depth. 



The optimal number of passes are obtained counting non-zero variables in the solution and 
associated depths are declared as optimal. The speeds and feeds corresponding to optimal 
depths are already available from single pass solutions. 

The common practice has been to use small depth of cut in finish pass and to remove 
the remaining depth of cut in a number of rough passes of equal depth yielding minimum 
production <.ost. I he multipass turning problem for all feasible total number of passes 
are formulated. T he minimum production cost corresponding to each such formulation is 
obtained and the lowest among these formulations is selected as optimal production cost 
and correspondiiag total number of passes are declared as optimal number of passes Shin 
and Joo [1992] have suggested that the small value for finish pass depth should be kept at 
the least allowable value of depth of cut to get the least value of production cost in finish 
pass. 

In this chapter five examples have been considered for the purpose of illustration. The 
first two exampkis show the use of integer programming model and compare the results ob- 
tained using iut;<^g(^r i>rogramming model and commonly practised strategy as suggested by 
Shin and Joo [1992]. T'h<; second example has also been solved using dynamic programming 
to compare with t.lu' i-('kuU,s obtained previously using integer programming and commonly 
practised strategy. 'I'ht'. saim; example is called example 3 when it is solved on the basis 
of varying workpiccx’; diann^ter using dynanoic programming model. These three examples 
have been solved at a given tool life in all passes in order to compare the results of three 
models viz. Shin and Joo’s model, integer programming model and dynamic programming 
model. The dynamic programming model as developed in the present work is able to handle 
varying workpiece diameter and non-uniform tool life in all passes. Example 4 has been 
taken as a general one and it has been found that, in general, the optimal tool life for finish 
pass is different from that obtained for rough passes. The fifth example has been taken to 
compare computer time for three other nonlinear methods viz. SUMT with DFP, GRG and 
SQP with dynamic i)rogramming model developed in the present work. 

4.2 Integer Programming Model 

This model is based on s<ime workpiece diameter in rough as well as finish passes. Further, 
it has been assumed that the optimal tool life for all passes are equal. The production costs ; 

for various depths of cut has been obtained on this assumption. The integer programming ; 

model works in two phases as described in section 3.4.1. In the first phase minimum 
production cost for a .sxu'ies of feasible depths of cut in finish as well as rough passes have been 
obtained. An integ<u’ {>rogramming problem equivalent to multipass turning is created in ^ 

the second phases and this is solved using LINDO software. The following two examples have | 

been solved using int<^g(u program ming model and the resulting depth of cut distribution ^ 



yields iiuiiiniuiM piodiK.tion cost in ixinjoritv of caso<^ TKa * f 

ha, be,,, solved „si.„„ L.NDO seftw ee wM b “ 

,K,.-fe.s s„l.dK.„s. such cases the integer programmrng fornrnlat™ is L able 

roach optnnal prodnot.on cost whereas dynamic programming formulation yields nnnime 

prodnction <^ost, 

In the proposed ...Odel, it has been assumed that the workpiece diameter is same in : 
the passi's and l.oldi rough and finish pass tools have the same tool life. These assumpti, 
have been ntatle to make problem formulation easy. The proposed model gives the value 
optimum aumlx'r of passes using branch and bound method. It has been found that t 
software LINDO sometimes fails to terminate at the optimal solution. Moreover the dep 
of cut values obta,iue(l wheu applied in different order will actually result in different coi 
when ,s(v>n iu eonjunetiou with the varying diameter model. The model is silent about t 
order in which t.lu'.se oi.timal <lepth of cut should be applied. 


4.2.1 Exainplu! 1 

Th(', multipa.s.s cylindrical turning example given by Shin and Joo [19921 is considered 
thi.s ('xaniph'. I he data its u.sed by Shin and Joo [1992] are given below. 

• IJlank (iiaiiK'tf'r 90 inm. 

• Lengtli of th<^ woi-ki)i('ce ~ 300 mm. 

• lotal depth of cut to he rt'moved = 6 mm. 

• Allowable cutting specui range, Vi min == 5 m/min; Uj = 500 m/min. 

• Allowable cutting feed range, fi min = 0.1 mm/rev; fi max = 0.9 mm/rev. 

• Allowable cutting depth range, di min = 1-0 nam; di max = 3.0 nun. 

• Constant for tool life equation, = 6 x lO'h 

• Ex[)onent,s for tool life ecpiation, ao,' = 5.0; /3o» = 1-75; 7oi = 0.75 or nu = 0 

nai = 0.35; = 0.15. 

• Allowal)le cutting tool life range, Ti min = 25 min; Tj max = 45 min. 

• Con.stant ftu- cutting force equation, K'n = 1.059 

Exponents for cutting force equation, On = 0, /3 h = 0.75; 7 i, = 0.95. 

• Maximum allowable cutting force, Fi max = 1.062 kN. | 

• Cutting j)ower e(|uation constant, K^i = 0.02076. j 


• Exponents for cutting power equation, asi = 1.0, p 2 i = 0.75; 73 ^ = 0.95. 

• Maximum allowable cutting power, P; mm = 5 kW. 

• Constant for surface finish equation, = 0.125/r; nose radius r = 1.2 mm 

• Exponents for surfac<' finish equation, = 0; = 2.0. 73 ^ = 0. 

• Maximum allowabki surface roughness for finish and rough passes Pq mm = 10 /rm; 

Ri max = 100 /i,m. 

• Constants for , tool approach and detraction time, = 7 x 10“^; h 2 i = 0.3 min. 

• Time components, = 0.75 min/piece; = 1.5 min/cutting edge. 

• Cost components, ku - $ 0.5/min; ku — $ 2.5/cutting edge. 

# 

For the above data (Xpi. (3.4) gives ^40^ = 0.249, assuming Ti = 25 min (test point), 
hitai = 0.255 and kiitp\ ~ 0.375 has to be added to the final solution once only. 

For total depth of cut valiu's upto 10.0 mm, the maximum number of rough passes are 
assumed to be 4, This is b(!cau.se the minimum depth of cut which can be removed in finish 
pass is 1 mm and therfdorc; mininmm number of rough passes required for Shin and Joo’s 
method are 3. Thc'-rofore the value of assumed maximum nmnber passes has been taken to 
be more than 3. Thus, integer programming problem can be formulated even for a larger 
number. The number of integer variables is more when a large value for maximum number 
of passes is selected. The restdt of integer programnnng problem remains unaffected with 
the choice of this number and therefore the smaller numbers should be preferred. 

The value of mi for i = 0, 1 n is taken as 20 for the generation of depth of cut series. 

Therefore, 

diQ — 1.0 mm 

dix = 1.0 + (3.0 ~ 2.0)i = 1.1 mm 
di 2 = 1.2 nun 
diz = 1.3 mm 


dm = 3.0 mm. 


The minimization of total cost Cto can be stated as 
minimize 

i=0 


( 4 . 1 ) 




(4.2) 


where the value of Cun for finish or rough pass as 


given by eqn. (3.3) is 


= + 0.255 

The cost foi finish pass and rough pass are minimized 


for that pas,, usin,; I h.' „».|,h.,.l„l„gy give,, ,eetion 3 2 

I j r 1 ’ J-ue two phases of the Droblpm 

can now l)e (uniiK'd. ^ px^oiem 


Phase 1 


The cost minimization for finish pass may be stated as 

minimize 


a,* = 0.249/„;«“,iS,>= + 0,255 

subject to th('. following constraints: 

(i) Cutting force constraint. 

(4.3) 

"•■'■"'C’C < 1 

(ii) Cutting pow(>r constraint, 

(4.4) 

0.4152x10 V,OC<l 

(4.5) 

(iii) Surface finish constraint 

/(), < 0.31 

(iv) Bounds on sixshI, feed and depth of cut 

, (4.6) 

5 < % < 500 

(4.7) 

0-1 < /o, < 0.9 

(4.8) 


Combining eqns. (4.6) and (4.8) the constraints for feed and depth of cut assume the form 


0.1 < /oi < 0.31 (4.9) 

1.0<^io,<3.0 ■ (4.10) 

The maximum ix'rmissihle feixi under the constraints to obtain the minimum production 
cost at the t,<!st {)oiat (25 min) is obtained using steps 1 to 3 described in section 3.2. There- 
3'fter the minimum production co.s’t and the speed corresponding to minimum production 
cost are ob<,ain<xi using stcj>s 4 and 5. In addition to this, the values for variable production 
cost per unit diamidcu- and variable production cost per unit diameter per unit 

length are also crahnilated for a given value of depth of cut do,. These values for all 

feasible depths of cut for finish pass are listed in Table 4.1, 


All constraints art' same for rongh passes except the value of 7? ..a .n 

max ana tne surface 

finish constraint takes the form 

fi, < 0,!I79 jjj 

The conetramt of aurfu,-,. (inieU therefore does not make any change on the bounds of feed 
and the production cost minimization problem for any rough pass is stated as 
minimizt^ 

C,„ = ().24<)/r»-“<i»,'= + 0.255 

subject to the following constraints: 

(i) Cutting force constraint 

< 1 

(ii) Cutting pow<'r constraint 

().4152 X 10 < 1 

(iii) Bounds on sp<'<>d, hs'd and dejith of cut 


(4.12) 

(4.13) 

(4.14) 


< -"it, < TiOO 
0.1 .< fi„ < 0.9 
l-O < < 3.0 


(4.15) 

(4.16) 

(4.17) 

The miiximum ix'rmi.ssibh' h'cd under the constraint for rough pass is obtained in the 
manner similar to that, lor linish pass. This feed value minimizes the TniniTmun production 
cost iy. In addition to the minimum production cost, the minimum variable production 
cost per unit diam(‘t,<'r ( tlie minimimi variable production cost per unit diameter 

per unit length and tlm cutting speed are obtained as for finish pass for each 

value of depth of cut. Tlu'st^ valutis for each value of depth of cut dig are listed in Table 4.2. 


Phase 2 


In this phas(' th<' al>ov<' probhun is formulated as a linear integer programme as described 
in phase 2 of sf-ctiou 3.4.1 and is .solved irsing LINDO software. The results obtained are 
given in Tahh'. 4.3. 'I'Ik' results for the same problem have also been obtained using Shin 
and Joo’s nuithod. d'hc.se are also given in Table 4.3. The integer programming method 
indicates that one rough pass and one finish pass are required. The depth of cut for both 
the passes is 3.0 mm and the minimum production cost is $ 1.94/piece. The method of Shin 
and Joo [1992] rcxpiirf's two rough passes of 2.5 mm and one finish pass of 1.0 mm depth 
and the corresponding minimum production cost is $ 2.385/piece. Therefore the integer 
programming method d(4;ermines lower production cost than that obtained by Shin and 
Joo. The total nuinbej- of pass<'s r<^quired in the proposed integer programming method to 
also lower or equal to that found using Shin and Joo’s method. 



Effect of Varying Total Depth of Cut on Minimum Cost 

of Shin and loo a.,. a, in Tablo 4.S for conapariaon poroses. For ar<ieXZZ 

to proposed .nodal roHnlts ... I„w,.r production coat than the method of Shin and Joo ^ 

„ptim.m. mnnl-cr o p..a,.a in the propo.«i method are also either lower or equal to tha 

found by Shux and Joo. " ° 

4.2.2 Example 2 

In example 1 , the depth of cut for both finiah and rough passes has been taken in thd range 
of 1.0 to ,i.O mm. However, li'oin the practical point of view this depth of cut range tor the 
anish pass see.n.s to be .... the higher .side and new range is taken as 0.3 to 1.2 mm (Hindnia 

, I , i)nun inin auu _ 1.2 mm. The depth of cut range for the 

rough pas.s aa<l ot.hor dat.a ar<> the .saim^. Th(^ two phases for the solution are described 
below 


Phase 1 


For da - 10.0 uuu, t,h(‘ maximum miml>or of passes io again taken as 4. The value of mo 
is taken as 8 for j-vm'ratiug t he <lt>i)th of cut series for finish pass, i.e., 
diH) “• O.J mm 
</()! 0.4 mm 


<k<) = 1-2 mm. 

The minimum pioductioix <-o,st for finish pass, (7^*^ q^, is found along with the optimum cut- 
ting speed feed for g = (), 1 — m.o. In addition to this the values for CjZr oq 
^umir Oil al.so culculatf'd cornssponding to each dog. These are listed in Table 4.4. 
The vahu's ol ■nt, ar<' t.akfui a.s 20 as in example 1 and the same dig series is gener- 
ated. Further, th<-, values (J,„ ami fl^ have not been calculated for this series. They 

are already availablf- in 'Fable 4.2. 

Phase 2 

The above problem is foi-mulat<'<l as an integer programming problem and is solved using 
LINDO softw;u-e as for t'xamph' 1. 'I'h<‘ results for this example using the integer program- 
Kdng method and Shiji and Joo’s Method are presented in Table 4.5. 



For da = 6.0 to 7.0 mm, the proposed integer programming model gives reduced number 
of passes and rcaluceci minimum production cost than those obtained using Shin and Joo’s 
method (or commonly pract,is(>d .strategy). For dty = 7.5 mm the proposed model does not 
achieve minimum production cost, probably because Lindo software sometimes terminates 
at a non-optimal and/or non-h-asibh' .sohition. It will be noted that the solutionis feasible 
but it may b<' non opt.imal (s(>(' (li.scusHiou.s of dynamic programming model). The solution 
of Shin and Joo’s nu'thod is available in feasible discrete space set. 

The propos(Hl mod('l giv('s slightly higher production cost of $ 2.876/per piece as com- 
pared to $ 2.7()()/pi('c(> obtaimul using the other method. 

The propos(Kl mode.l giw^s slightly higher minimum production cost of $2.946/piece at 
dti = 8.0 mm. This tinw' th<' h'asible discrete space in the proposed model does not contain 
the minimum cost point, obtaimsl by the other method. The proposed model success- 
fully achi<w('s <'ith<'r ('(jual or naluced minimum production cost of $ 3.006/piece and $ 
3.020/pi('c<^, r('sp<'ct.ivf'ly for da 8.5 mm and da = 9.0 mm. Further, at dti = 9.5 mm 
or 10.0 mm t.lu' propos«><l model giws re<hiced number of passes and reduced mim'TnnTY^ 
production cost, (in <lisrr<>t.e feasil)h^ S}>ace) than those obtained using the other method. 
Thus, out of 9 cases consid<'r«Hl abov<', the. proposed model clearly gives minimum produc- 
tion cost at maximum of finishing depth range, i.e., 1.2 mm for 5 cases by reducing number 
of passes. For da 7.5 to 9.5 mm, the proposed model gives almost equal or reduced 
minimum pnxluct.ion cost for two ca.s(is and failed in remaining two cases. These two cases 
will be exarniiu'd in det.ail using dynamic programming model. 

4.3 Dynamic Programming Model 

The integer programming moch'l is ba.sed on constant initial workpiece diameter whereas 
the dynamic progrannning moikd de.scribed in section 3.4.2 is able to acconomodate varying 
workpiece diam(4,(us at (uich pfiss. It can also be adapted to constant initial workpiece 
diameter. Th<i solut.ion proexsinn^ for dynamic programming model is given below. 

4.3.1 Example 3 

In this (example the daf.a is f,hc sanw* a.s in example 2. It has been assumed that the tool life 
for both rough as wcdl as finish pass is the sjxme and is equal to 25 min. This assumption 
has been rnah; in ord<‘r f.o compare the results obtained using dynamic programming model 
with thos('. obt.aiiu'd pixwiously using integer programming model or Shin and Joo’s method 
where tool lift; in all passers has btxm assumed to be the same. 



Solution Procedure 


The initial workpiece diameter is 50 mm. If the total depth of cut to be removed is 6 mm, 
the final workpiece diameter will be 50— 2 x 6 = 38 mm. Thus 6 mm depth is divided into 60 
divisions of 0.1 mm (the number of divisions depend upon the desired accuracy). Similarly, 
the depth of cut ranges for finish as well as rough passes are converted into number of 
divisions as 3 to 12 and 10 to 30, respectively. 

Let a matrix element cost^q be defined as minimum production cost for machining total k 
divisions when the depth of cut corresponding to division is machined in rough pass. The 
value of q has been defined in section 4.2.1. The minimum cost for machining k divisions, 
costmin k is given as 

costmin 1 = minimum costkg for all q (4.18) 

costmin 1 and costmin 2 are given a very high value of 9999999 because finish pass of 0.1 and 
0.2 mm depths of cut are not feasible. 
costmin 3 to costmin 12 are Calculated as 

costmin k = (-02/; + 0 ( 4 - 3 ) + 0.255 (4.19) 

where the value of 0 ( 4 - 3 ) i® obtained from Table 4.1. Using the relationship 

costkq = (0.02k + Df)Ct^^ ig + costminii'>T^x) +0.255 (4.20) 

costkq for k = 13 to 60 and 5 = 0 to 20 are obtained provided q = (k -(3 + 10)) is non- 
negative. Here iq i® obtained from Table 4.2, index is given as 

index = k — (0.3 -f- q) (4.21) 

and costkq = 9999999 if term [A; - (3 -|- 10)] is negative, costmin k is then obtained and the 
value of q for which costmin k corresponds is called keyk- This key is used for backtracking 
the depth of cut distribution, costmin eo is foimd and keyeo gives the value of depth of cut 
for the first rough pass and 60 - keyeo gives the remaining depth of cut and keyeo-ksyeo ’^iii 
give the depth of cut for the second pass and so on. 

Observations 

The values for the minimum production cost for various depths of cut are given along with 
the depth of cut distribution for finish pass as well as rough pass in Table 4.5. The dynamic 
programming model has been adapted for constant diameter model using Dq in place of 
(0.02A; + Df) in eqns (4.19) and (4.20). The results for various depths of cut have been 
presented in Table 4.6 and it is clear that for all depths of cut ranging from 6 to 10 mm the 
minimum production cost obtained using dynamic progra mmi ng model is lower than that 
obtained using Shin and joo’s method or integer programming model. 
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It is important to note that for a depth of cut of 7.5 mm the integer progr amming model 
is not able to achieve lower minimum production cost than $ 2.872/piece as obtained by 
Shin and Joo’s method (see Table 4.5) but the dynamic progra mming model gives lower 
minimum production cost of $ 2.848/piece (see Table 4.6). Similarly, for a depth of cut of 8 
mm also when integer programming model is not able to achieve the minimum production 
cost, the dynamic programming model yields minimum production cost. The most probable 
reason for these kind of incidents is that LINDO software used in the second phase of integer 
programming model sometimes gives non-optimal results. It is also clear from the results 
given in Table 4.6 that the minimum production cost with varying workpiece diameter gives 
lower value than that obtained using initial workpiece diameter. Thus it is evident that 
dynamic programming model gives more accurate results than integer programming model. 
Further, integer programming model does not give any information about the order in which 
the depths of cut should be applied but the dynamic programming model gives the order 
in which these optimal depths of cut should be applied. 

4.3.2 Example 4 

All the three examples presented earlier were minimizing the production cost at tool life of 
25 min which remain the same for both finish and rough passes. The assmnption of equal 
tool life in all passes is, in general, not true but it has been used in these three examples to 
compare the results of dynamic programming model, integer programming model and Shin 
and Joo’s method. Therefore in this example the tool life and constraint data has been 
taken partially from Ermer and Kromodihardjo [1981] and the cost and time coefficients 
are more or less the same as in example 1. The example data are as follows: 


• Blank diameter = 50 mm. 

• Length of the workpiece = 300 mm. 

• Total depth of cut to be removed = 6 mm. 

• Allowable cutting speed range, Vi min = 5 m/min; Vi max = 500 m/min. 

• Allowable cutting feed range, fi min = 0.1 mm/rev; fi mxix — 0.9 mm /rev. 

• Allowable cutting depth range for finish pass, do min = 0.3 nun; do max = 1-2 mm. 

• Allowable cutting depth range for rough passes, di min = 1-0 mm; di max = 3-0 mm. 

• Constant for tool life equation, Koi = 1.40673 x 10®. 

• Exponents for tool life equation, aoi = 4.0; j3oi = 1.16; 7oi = 0.15. 
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• Allowable cutting tool life range, T; mm = 5 rnin; 'i- max = *)U min. 

• Cutting power equation constant, K'y = 0.499. 

• Exponents for cutting power equation, an = 0.95, /3n = 0.78; 7ii = 0.75. 

• Maximum allowable cutting power, Pi max = 20 kW. 

• Constant for surface finish equation, K^i = 0. 125 /C' .nose radius r = 1.2 mm. 

• Exponents for surface finish equation, a2i-d-5'l/32i = [•oo 4 - 72 i = 0.2..5 

• Maximum allowable surface roughness for finish and rough passes i?o max = 5 pm- 

Pi max ~ 25 /im. 

> 

• Constants for tool approach and detraction time, hu = 7 x 10"'*; /i 2 t = 0.3 min. 

• Time components, fpi = 0.75 min/piece; = 1-5 min/cutting edge. 

• Cost components, ku = $ 0.5/min; kn = $ 2.5/cutting edge. 

In this example the tool life, cutting power and surface finish equations have been taken from 
Ermer and Rrmodihardj o [1981], 

Observations 

The minimum production cost for different depths of cut ranging from 0.3 to 1.2 mm are 
listed in Table 4.7. The optimal tool life corresponding to various depths of cut changes 
from 15.78 to 38.85 min in finish pass depending upon the depth of cut used. 

The minimum production cost for various depths of cut ranging from 1.0 to 3.0 mm ir 
rough pass are given in Table 4.8. The optimal tool life corresponding to various depth; 
of cut are the same (59.97 min). It may be noted that this value is significantly differen 
than that obtained in the finish pass (15.78 to 38.85 min). Therefore the assumption o 
equal tool life in all passes made to simplify the solution of multipass turning problem wil 
lead to incorrect results. The minimum production cost in this case is $ 4.4457/piece base 
on variable workpiece diameter and without any restriction on tool life in any pass. Th 
optimum machining conditions are as following; 

= 1.2 mm; = 143.9 m/min; = 0.853 mm/rev; 

= 2.3 nam; = 83.93 m/min; = 0.231 mm/rev 

dT = 2.5 mm; = 83.87 m/min; = 0.214 mm/rev 



4.4 Comparison of Various Nonlinear Techniques 

The various nonlinear programming techniques for minimization of production cost have 
been described in section 3.5. These are named as sequential unconstrained ruiniinization 
technique (SUMT) with Davidon, Fletcher and Powell’s (DFP) method, generalized reduced 
gradient (GRG) method and sequential quadratic programming (SQP) nx't.hod. 

4.4.1 Example 5 

The objective function chosen by Chua et al. [1993] is minimum production time and oui} 
rough turning is considered. The production cost function results in production time func 
tion if overliead cost is taken as unity and machining as well iris tool edgt^ cost are takei 
to be zero. Further, the surface finish requirement is slackened or a Ugh value of Rr,ua 1 
selected for both rough and finish passes. The following data are ust'd: 

The data used for machining of Rochling T4(C 0.45%, Si 0.25%, Mil 0.70%) with SANDVIU 42 
TiN-coated carbide insert (SNUN 120408) by Chua et al. [1993] arc given below. 

• Blank diameter = 200 mm. 

• Length of the workpiece = 400 mm. 

• Total depth of cut to be removed = 4 mm. 

• Allowable cutting speed range, U{ min — 150 m/min; Vi „mx — 212 m/min. 

• Allowable cutting feed range, fi min = 0.2 mm/rev; /; ~ 0,35 rnm/rov. 

• Allowable cutting depth range, di min - 1-0 mm; di max = 2.0 rnm. 

• Constant for tool life equation, Koi = 8.81 x 10^. 

• Exponents for tool life equation, ctoi = 2.894; poi = 0.493; 7oi = 0. 

• Allowable cutting tool life range, Ti mir. = 27.373755 min; Ti = 98.165 min. 

• Constant for cutting force equation, = 10.1 

• Exponents for cutting force equation, an = -.0.171, /),, = 0.641; 7,, ^ 0.852. 

• Maximum allowable cutting force, Fi max - 12 kN. 

• Cutting power equation constant, K^i = 0.179. 

• Exponents for cutting power equation, 0:2.' = 0.814, p^i 0.634; ().845. 

• Maximum allowable cutting power, max = 3 kW. 



• Constant for surface finish equation, — 0.125/r; nose radius r = 1.2 mm. 

• Exponents for surface finish equation, 0:31 = 0; /?3i = 2.0. qsi = 0. 

• Maximum allowable surface roughness for finish and rough passes, Rq max = 2500 /xm 
Rimax = 2500 nm. 

• Constants for tool approach and detraction time, hu = 0.0; /i2i = 0.0 min. 

• Time components, tpi = 1.0 min/piece; td = 1-0 min/ cutting edge. 

• Cost components, Jsii = $ 1.0/min; ka = $ 0.0/cutting edge. 

(It may be noted that this is a production time minimization problem because Ku and 
are unity and ku are zero). 

4.4.2 Observations 

All the above mentioned techniques require the number of passes to be known for the for- 
mulation of optimization problem. Umin = 2 Umax = 4 for the above problem. Therefore the 
multipass turning model axe to be formulated and solved for 2, 3 and 4 number of passes. 
Otherwise three problems of similar nature have to be solved as suggested by Chua et al. 
[1993]. In the first problem, the number of passes are treated as continuous variable and an 
optimum continuous approximation of the mmober of passes has been obtained as 2.58. The 
next two problems are formulated for two nearest integer to the continuous optimum num- 
ber of passes, i.e., 2 and 3. The first two methods require feasible starting point for which a 
separate optimization problem is to be solved. SQP may be started from an arbitrary point 
and dynamic programming alongwith single pass solution methodology does not require any 
starting point. SUMT with DFP method is dependent on optimization parameters. Sim- 
ilarly, GRG is also dependent on optimization parameters and convergence of the method 
depends upon the initial partition of basic and non-basic variables. On the otherhand, no 
such parameter are required for SQP method. The dynamic programming method neither 
requires any optimization parameter nor it requires any knowledge of number of passes. 
Rather the dynamic programming method determines the number of passes uniquely. The 
optimum distribution of depth of cut is also successfully determined. Therefore this method 
should be preferred over the other three because of three reasons: (i) the number of passes 
are determined uniquely, (ii) no prior knowledge of feasible starting point is required and 
(iii) no optimization parameter is required. The computation time which is the sum of user 
time and system time have been compared for the three non-linear techniques alongwith 
the dynamic programming model using single pass solution methodology. It has been found 
that SQP requires the least time, that is approximately 3 x 0.1 = 0.3 seconds. This method- 
ology may not be available to all manufactures. The dynamic programming methodology 
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requires 0.9 seconds but is easy to implement. GRG method ranks the third and it requires 
3 X 0.4 = 1.2 seconds and SUMT with DFP requires 3 x 0.8 = 2.4 seconds. It may be noted 
that the last two methods GRG and SUMT with DFP require initial feasible solution also. 
The dynamic programming method should, therefore, be preferred over other methods. 
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S.No. 

doq 

mm 

^0? 

m/min 

foq 

mm/rev 

^tovar Oq 

$/ piece 
10-3 

^tovar Oq 

$/ piece 

10-3 

^toOq 

$ /piece 

0 

1.0 

179.5 

0.31 

10.66 

3.553 

0.788 

1 

1.1 

177.0 

0.31 

10.82 

3.607 

0.796 

2 

1.2 

174.7 

0.31 

10.96 

3.653 

0.803 

3 

1.3 

172.6 

0.31 

11.08 

3.693 

0.809 

4 

1.4 

170.7 

0.31 

11.22 

3.740 

0.816 

5 

1.5 

168.9 

0.31 

11.34 

3.780 

0.822 

6 

1.6 

167.3 

0.31 

11.44 

3.813 

0.827 

7 

1.7 

165.8 

0.31 

11.54 

3.847 

0.832 

8 

1.8 

164.4 

0.31 

11.64 

3.880 

0.837 

9 

1.9 

163.0 

0.31 

11.74 

3.913 

0.842 

10 

2.0 

161.8 

0.31 

11.84 

3.947 

0.847 

11 

2.1 

160.6 

0.31 

12.00 

4.000 

0.855 

12 

2.2 

159.5 

0.31 

12.00 

4.000 

0.855 

13 

2.3 

158.4 

0.31 

12.08 

4.027 

0.859 

14 

2.4 

157.4 

0.31 

12.16 

4.053 

0.863 

15 

2.5 

156.4 

0.31 

12.24 

4.080 

0.867 

16 

2.6 

155.5 

0.31 

12.30 

4.100 

0.870 

17 

2.7 

154.7 

0.31 

12.38 

4.127 

0.874 

18 

2.8 

153.8 

0.31 

12.44 

4.147 

0.877 

19 

2.9 

153.0 

0.31 

12.50 

4.167 

0.880 

20 

3.0 

152.2 

0.31 

12.58 

4.193 

0.884 


Table 4.1; Optimal cutting conditions and costs when the total stock is removed in a single 
finish pass at Tq = 25 min (Example 1). 
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S.No. 

diq 

mm 

m/min 

f* 

J iq 

mm /rev 

^tavar iq 

$ 10“^ /piece 

^tovar iq 

$ 10“®/piece 

^toiq 

$/piece 

0 

1.0 

123.6 

0.9 

5.34 

1.780 

0.522 

1 

1.1 

121.8 

0.9 

5.40 

1.800 

0.525 

2 

1.2 

120.3 

0.9 

5.48 

1.827 

0.529 

3 

1.3 

118.8 

0.9 

5.54 

1.847 

0.532 

4 

1.4 

117.5 

0.9 

5.60 

1.867 

0.535 

5 

1.5 

116.3 

0.9 

5.67 

1.887 

0.538 

6 

1.6 

115.2 

0.9 

5.72 

1.901 

0.541 

7 

1.7 

114.1 

0.9 

5.78 

1.927 

0.544 

8 

1.8 

113.2 

0.9 

5.82 

1.940 

0.546 

9 

1.9 

112.3 

0.9 

5.88 

1.960 ■ 

0.549 

10 

2.0 

111.4 

0.9 

5.92 

1.970 

0.551 

11 

2.1 

111.1 

0.889 

6.00 

2.000 

0.555 

12 

2.2 

112.6 

0.838 

6.28 

2.093 

0.569 

13 

2.3 

114.1 

0.792 

6.56 

2.187 

0.583 

14 

2.4 

115.5 

0.750 

6.84 

2.280 

0.597 

15 

2.5 

116.9 

0.712 

7.12 

2.370 

0.611 

16 

2.6 

118.3 

0.678 

7.40 

2.407 

0.625 

17 

2.7 

119.6 

0.646 

7.68 

2.560 

0.639 

18 

2.8 

120.9 

0.617 

7.96 

2.650 

0.653 

19 

2.9 

122.2 

0.590 

8.24 

2.747 

0.667 

20 

3.0 

123.4 

0.565 

8.52 

2.840 

0.681 


Table 4.2: Optimal cutting conditions and costs when the total stock is removed in a single 
rough pass at Tj = 25 min (Example 1). 
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Proposed model 

Shin and Joo’s model 

S.No 

dti 


df 

cif 

df 

n°pt 

QOpt 

df 

df 


QOpt 


mm 

mm 

mm 

mm 

mm 


$ /piece 

mm 

mm 


$/piece 

1 

6.0 

3.0 

3.0 



1 

1.940 

1.0 

2.50 

2 

2.385 

2 

8.0 

3.0 

2.1 

2.9 


2 

2.481 

1.0 

2.33 

3 

2.927 

3 

8.5 

3.0 

2.5 

3.0 


2 

2.551 

1.0 

2.50 

3 

2.996 

4 

9.0 

3.0 

3.0 

3.0 


2 

2.611 

1.0 

2.67 

3 

3.068 

5 

9.5 

2.9 

1.9 

2.8 

1.9 

3 

3.005 

1.0 

2.83 

3 

3.134 

6 

10.0 

3.0 

2.1 

2.8 

2.1 

3 

3.022 

1.0 

3.00 

3 

3.206 

7 

10.5 

2.7 

2.0 

2.9 

2.9 

3 

3.134 

1.0 

2.50 

4 

3.539 

8 

11.0 

3.0 

2.9 

2.1 

3.0 

3 

3.162 

1.0 

2.50 

4 

3.607 


Table 4.3: Comparison of results obtained for the proposed integer programming model and 
Shin and Joo’s model at Tj = 25 min (Example 1). 


S.No. 

doq 

mm 

^0? 

m/min 

foq 

mm/rev 

^tovar Oq 

$ /piece 

^tovar 0^ 

$ /piece 

n* 

^tovar Qq 

$ /piece 

0 

0.3 

215.0 

0.31 

8.9 

2.967 

0.700 

1 

0.4 

205.9 

0.31 

9.3 

3.100 

0.720 

2 

0.5 

199.2 

0.31 

9.6 

3.200 

0.735 

3 

0.6 

193.8 

0.31 

9.88 

3.293 

0.749 

4 

0.7 

189.4 

0.31 , 

10.10 

3.367 

0.760 

5 

0.8 

185.6 

0.31 

10.32 

3.440 

0.771 

6 

0.9 

182.4 

0.31 

10.50 

3.500 

0.780 

7 

1.0 

179.5 

0.31 

10.66 

3.553 

0.788 

8 

1.1 

177.0 

0.31 

10.82 

3.607 

0.796 

9 

1.2 

174.7 

0.31 

10.96 

3.653 

0.803 


Table 4.4: Optimal cutting conditions and costs for modified finish pass range at To — 25 
min (Example 2). 
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Proposed model 

Shin and Joo’s model 

S.No 

dti 



df 

df 


QOpt 


df 


QOpt 


mm 

mm 

mm 

mm 

mm 


$/piece 

TYITTI 

mm 


$/piece 

1 

6.0 

1.2 

2.7 

2.1 


2 

2.372 

0.3 

1.90 

3 

2.722 

2 

6.5 

1.2 

3.0 

2.3 


2 

2.442 

0.3 

2.07 

3 

2.736 

3 

7.0 

1.2 

3.0 

2.8 


2 

2.512 

0.3 

2.23 

3 

2.796 

4 

7.5 

0.3 

2.7 

2.5 

2.0 

3 

2.876 

0.3 

2.40 

3 

2.866 

5 

8.0 

0.3 

2.7 

3.0 

2.0 

3 

2.946 

0.3 

2.57 

3 

2.936 

6 

8.5 

0.3 

2.7 

3.0 

2.5 

3 

3.006 

0.3 

2.73 

3 

3.020 

7 

9.0 

0.3 

3.0 

2.7 

3.0 

3 

3.076 

0.3 

2.90 

3 

3.076 

8 

9.5 

1.2 

3.0 

3.0 

2.3 

3 

3.123 

0.3 

2.30 

4 

3.407 

9 

10.0 

1.2 

3.0 

3.0 

2.8 

3 

3.193 

0.3 

2.43 

4 

3.477 


Table 4.5: Coraparison of results obtained for the proposed integer programming model and 
Shin and Joo’s model for modified jfinished pass range at 71 = 25 min (Example 2). 





DP model var 

dia 


DP model constant dia 

S.No 

dt 

mm 

df 

mm 

df 

mm 

df 

mm 

df 

mm 


QOpt 

$ /piece 

df 

mm 

df 

mm 

df 

mm 

mm 

n°pt 

QOpt 

$ /piece 

1 

6.0 

1.2 

2.4 

2.4 


2 

2.234 

1.2 

2.1 

2.7 


2 

2.372 

2 

6.5 

1.2 

2.7 

2.6 


2 

2.286 

1.2 

3.0 

2.3 


2 

2.442 

3 

7.0 

1.2 

2.9 

2.9 


2 

2.337 

1.2 

3.0 

2.8 


2 

2.512 

4 

7.5 

1.2 

2.1 

2.1 

2.1 

3 

2.629 

1.2 

2.1 

2.1 

2.1 

3 

2.843 

5 

8.0 

1.2 

2.3 

2.3 

2.3 

3 

2.676 

1.2 

2.6 

2.1 

2.1 

3 

2.913 

6 

8.5 

1.2 

2.5 

2.4 

2.4 

3 

2.722 

1.2 

3.0 

2.2 

2.1 

3 

3.020 

7 

9.0 

1.2 

2.6 

2.6 

2.6 

3 

3.676 

1.2 

3.0 

2.7 

2.1 

3 

3.053 

8 

9.5 

1.2 

2.8 

2.8 

2.7 

3 

2.810 

1.2 

3.0 

3.0 

2.3 

3 

3.123 

9 

10.0 

1.2 

3.0 

2.9 

2.9 

3 

2.853 

1.2 

3.0 

3.0 

2.8 

3 

3.193 


Table 4.6: Comparison of results obtained for dynamic programming model using variable 
and constant diameters for modified finish pass range at Tj = 25 min (Example 3). 
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S.No. 

d(jq 

mm 

m/min 

foq 

mm/rev 

n-i* 

J-Oq 

min 

^tovar Oq 

$ /piece 
10-2 

^taoar Oq 

$/piece 

10-2 

^tdOq 

$/piece 

0 

0.3 

143.9 

0.853 

15.78 

5.437 

1.813 

0.515 

1 

0.4 

133.04 

0.709 

19.01 

6.702 

2.234 

0.517 

2 

0.5 

174.7 

0.614 

21.98 

7.920 

2.640 

0.518 

3 

0.6 

172.6 

0.546 

24.74 

9.102 

3.034 

0.519 

4 

0.7 

170.7 

0.494 

27.36 

10.026 

3.419 

0.520 

5 

0.8 

168.9 

0.454 

29.83 

11.388 

3.796 

0.521 

6 

0.9 

167.3 

0.421 

32.21 

12.500 

4.167 

0.523 

7 

1.0 

165.8 

0.393 

34.49 

13.598 

4.533 

0.524 

8 

1.1 

164.4 

0.370 

36.68 

14.680 

4.894 

0.525 

9 

1.2 

163.0 

0.350 

38.85 

15.750 

5.252 

0.526 


Table 4.7: Optimal cutting conditions and costs when the total stock is removed in a single 
finish pass (Example 4). 
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S.No. 


jr* nnt 

Jiq 

mm/rev min 


84.60 

84.50 

84.45 

84.39 

84.33 

84.27 

84.22 

84.16 

84.13 

84.08 

84.04 

84.01 

83.97 

83.93 

83.89 

83.87 

83.84 

83.81 

83.77 

83.75 

83.72 


0.51 

0.466 

0.429 

0.397 

0.370 

0.347 

0.326 

0.308 

0.292 

0.277 

0.264 

0.252 

0.241 

0.231 

0.222 

0.214 

0.206 

0.199 

0.192 

0.185 

0.180 


>7 
)7 
17 
17 
•7 
•7 
i7 
'7 
7 
7 
7 
7 
7 
7 
7 
7 

59.97 

59.97 

59.97 
7 
7 


C*' 

^tcrvar iq 

$ /piece 
10-3 


12.105 

13.264 

14.419 

15.570 

16.717 

17.861 

19.00 

20.141 

21.28 

22.401 

23.540 

24.670 

25.790 

26.917 

28.039 

29.159 

30.277 

31.394 

32.508 

33.622 

34.73 


/-I*" 

^tcwar iq 

$ /piece 
10-3 


4.035 

4.421 

4.807 

5.190 

5.572 

5.954 

6.334 

6.714 

7.092 
7.470 
7.846 
8.222 
8.598 
8.972 
9.346 
9.720 

10.092 
10.465 
10.836 
11.207 
11.578 


$/piece 


Table 4.8: Optimal cutting conditions and costs when the total stock is removed in a single 
rough pass (Example 4). 







S.No. 

Characterstics 

SUMT 

GRG 

SQP 

DP using 



DFP 



single pass 

1 

Min. prod, time 

10.467 

10.391 

10.37 


2 

Real time 

1.3 

■||||Q[|| 

0.2 

1.1 


User time 

0.8 


0.1 

0.9 



0.0 


0.0 

0.0 

3 

Feasible initial 

Yes 

Yes 

No 

No 


point required 
Initial speed 

170.0 


170.0 




170.0 


170.0 

~ 



170.0 


170.0 

- 


Feed 


0.25 

0.25 





0.25 

0.25 

- 




0.25 

0.25 

~ 


Depth of cut 

1.33 

1.33 

1.33 

- 



1.33 

1.33 

1.33 

- 



1.34 

1.34 

1.34 

1 

4 

Optimal speed 

212.04 

212.0 

212.0 




211.23 

212.0 

212.0 




211.28 

212.0 

212.0 



Feed 

0.3487 

0.3499 

0.3499 

0.2512 



.3480 

0.3499 

0.3499 

0.2512 



.3486 

0.3499 

0.3499 



Depth of cut 

1.449 

1.3455 

1.55 

2.0 



1.3499 

1.3248 

1.45 

2.0 



1.1977 

1.3297 

1.00 


5 

No. of passes 
required for 
problem 
formulation 

Yes 

Yes 

Yes 

No 

6 

Passes taken to 

3 

3 

3 

Not 


assess comp, time 




required 


Optimal passes 
obtained directly 

No 

No 

No 

Yes (2) 


Procedure to 

All 

All 

Treating no. of pass 



obtain optimal 

possible 

possible 

continuous variable 

Uniquely 


number of passes 

no. of 

no. of 

and solving two more 

obtained 



passes 

passes 

problems arround cont- 




(2,3,4) 

(2,3,4) 

inuoxis no. of passes 






2.58, (2,3) 



Table 4.9; Optimal cutting conditions, and costs and times for various nonlinear techniques 
and dynamic programming method (Example 5). 
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Chapter 5 

Turning of Stepped Shafts 


5.1 Introduction 

In the previous chapters the emphasis has been on optimization of single as well as multipass 
straight turning operation. It has been found in section 4.3 that DP approach does not 
require the knowledge of optimal number of passes in order to make multipass turning model. 
Instead the optimal number of passes are uniquely determined using DP approach which 
requires repetitive single pass solutions. The single pass solution methodology discussed 
in section 3.2 has advantage over other techniques that the optimal machining conditions 
are determined without the knowledge of feasible/nonfeasible starting point and without 
experimentation with optimization parameters. Thus, DP alongwith the single pass solution 
methodology discussed in the present work has its own advantages and has been selected to 
find the optimal clubbing of machinable volumes in vertical, horizontal or circular manner 
for stepped shafts. Figure 5.1 shows an example of rough turning of a stepped shaft. 
The various machinable volume segments (1), (2) and (3) have been obtained by drawing 
horizontal and vertical lines from the ends of cylindrical volrunes. These volumes are defined 
as machinable volume, seed volume, and compound volumes. 

Machinable Volume 

The smallest volrnne segment enclosed by two consecutive horizontal and two consecutive 
vertical lines is called the machinable volume. In Fig. 5.1, (1), (2) and (3) are machinable 
volumes. 

Seed Volume 

It is a machinable volume which has at least one surface as machined boundary of shaft. 
Volumes (1) and (5) are seed voliunes. 

Compound Volume 

This is a volume consisting of a seed and a few non-seed volumes, which can be machined 
together. For example volumes (1,3) and volumes (2,5) are compound vplumes. In order 
to make feasible sets of compound voltunes a seed volume is selected. Starting with the 
seed volume, say (5), first vertical clubbing is performed and volumes (5,6) and (5,6,7) 
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are generated. Thereafter, horizontal clubbing of volumes (5,2) is performed. The circular 

clubbing is generated from small to larger compound volumes (5, 2, 3,6) and (5,2,3,4,6,7). 

It is interesting to n^ that grouping of volumes (1,2) and (5) volumes (5,2) and (1) are 

two alternative/ tor macmmng. These alternative . will give different minimum production 
^ ^ 

costs. Therefore it becomes an interesting decision to select optimal grouping of machinable 
volumes so that the minimum production cost is obtained. For all seed and compound vol- 
umes, the optimal number of passes has to be determined for optimal grouping of machinable 
vol u mes to machine the shaft. 

The horizontal grouping of machinable volumes has been considered only in recent lit- 
erature for turning of rough stepped shafts [Jha, 1996; Prasad et al., 1997]. In a similar 
work for the milling process, the machinable volumes have been grouped using a heuristic 
based on ascending machining cost per unit volmne [Yellowley and Fisher, 1994]. A linear 
integer model has been suggested in this chapter to find the miuimum production cost and 
optimal grouping of machinable volumes for stepped shafts. It has been verified that the 
suggested model gives better mi nimum production cost than that obtained using heuristic 
method or common practice. 

5.2 Optimal Grouping of Machinable Volumes 

In order to make optimal grouping of machinable volinnes, all seed volumes are recognized 
and compound volinnes corresponding to all seed volumes are generated. These seed and 
compound volumes are candidates to be considered for grouping. Let a candidate volume be 
represented as Vg and the associated minimnni production cost obtained using DP alongwith 
single pass methodology be Ctog- The candidate volume is either a seed or a compound 
volume. A matrix A is defined to represent the composition of machi n ing volumes such 
that a^g is one if machinable volume e belongs to candidate volume g and is zero otherwise. 
The minimization of production cost for turning of a stepped shaft is defined as 
minimize 

j:c^^g (5.1) 

9 

where xg is a binary variable. Its value is one if g^^ volume is selected in optimal grouping 
and is zero otherwise, subject to 

J2^e,xg>l for all e (5.2) 

9 

The above integer programming formulation is solved using LINDO optimizer. 
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5.2.1 Example 1 

A stepped shaft shown in Fig. 5.1. has six cylindrical faces. The manufacturing of this 
shaft has been performed with a blank. It may be noted that surfaces (i), (ii) and (iii) 
can be machined while tool approaches from the right and surfaces (v) and (vi) can be 
machined while tool approaches from the left side. Surface (iv) can be machined from tool 
approaching from either side. 



Figure 5.1: Rough turning of a stepped shaft from 40 mm diameter blank 

In turning of stepped shafts, various small volume segments can be formed by drawing 
vertical lines from the ends of different cylindrical surfaces and horizontal lines from ends 
of various facing surfaces. The various volume segments thus formed are denoted from 1 to 
15. The following are the example data: 

• Blank diameter = 40 mm. 

• Allowable cutting speed range, Vi min = 5 m/min; Vi max = 500 m/min. 

• Allowable cutting feed range, fi min = 0-1 mm/ rev; / ; max = 0.9 mm/rev. 

• Allowable cutting depth range, dj min = 1-0 mm; dj max = 3.0 mm. 
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• Constant for tool life equation, = 6 x 10^^. 

• Exponents for tool life equation, aoi = 5.0; /3oi = 1.75; 7ot = 0.75. 

• Allowable cutting tool life range, Tj rnin = 25 min; Ti ^ = 45 min. 

• Constant for cutting force equation, = 1.059 

• Exponents for cutting force equation, au = 0, = 0.75; qn = 0.95. 

• Maximum allowable cutting force, Fi = 1.962 kN. 

• Cutting power equation constant, = 0.02076. 

• Exponents for cutting power equation, Q' 2 i = 1.0, ^ 2 i — 0.75; 72 i = 0.95. 

• Maximum allowable cutting power. Pi max = 5 kW. 

• Constant for surface finish equation, K'^i = 0.125/r; nose radius r = 1.2 mm. 

• Exponents for surface finish equation, ^zi = 2.0. 

• Maximum allowable smface roughness, Ri max = 100 /im. 

• Constants for tool approach and detraction time, hn = 7 x 10~^; hzi = 0.3 min. 

• Time components, tpi = 0.75 min/piece; td = 1.5 min/cutting edge. 

• Cost components, ku = $ 0.5/min; ka = $ 2.5/cutting edge. 

For the data, the constants calculated using eqn. (3.4) are: 

AOii = 0.249, assuming Tj = 25 min (test point) 
kiitpi = 0.375 


Solution Procedure 

There are six seed volumes in this example. These are volmnes (1), (5), (8), (10), (11) 
and (13). The compound volumes corresponding to each of these seed volumes are given in 
Table 5.1. The initial and final diameters and lengths of steps, mi n i mum production cost, 
volume and minimum production cost / volume are also given in the same table. There are 
15 machinable voltimes and 35 candidate volumes. The minimum production cost for each 
candidate volume have been obtained using DP alongwith single pass solution methodology 
and are given in Table 5.1. The integer progr ammi ng model formulated to obtain the 
minimum production cost as well as optimal grouping of machinable volumes is 
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nunimize 

0.3356x1 + 0.5154x2 + 0.6913x3 + 0.6999x4 + 0.3546x5 + 0.3772x6 + 0.5479x7 -f 0.3940x8 + 
0.4330x9 + 0.6185x10 + 0.2177x11 + 0.4172x12 + 0.2537x13 + 0.2797x14 + 0.4796x15 + 
0.5246x16 + 0.1580x17 + 0.2153x18 + 0.2458x19 + 0.2678x20 + 0.1986x21 + 0.2206x22 + 
0.2559x23 + 0.2589x24 + 0.2809x25 + 0.2779x26 + 0.3084x27 + 0.3305x28 + 0.3914x29 + 
0.4072x30 + 0.4408x31 + 0.4652x32 + 0.3417x33 + 0.6911x34 + 0.6997x35 
subject to 

xH- x2 + x3 4- x4 >= 1 

x2 + x3 + x4 + x8 + x9 + xlO >= 1 

x3 + x4 + x9 + xlO + xl4 + xl6 >= 1 

x4 + xlO + xl6 + x20 4- x25 4- x28 >= 1 

x5 + x6 4- x7 + x9 4- xlO >= 1 

x6 4- x7 4- x9 4- xlO 4- xl3 4- xl4 + xl5 4- xl6 >= 1 

x7 + xlO + xl5 H- xl6 + xl9 4- x20 4- x24 4- x25 + x27 + x28 >= 1 

xll 4- xl2 4- xl3 4- xl4 4- xl5 4- xl6 >= 1 

xl2 4- xl5 + xl6 4- xl8 4- xl9 + x20 4- x23 4- x24 4- x25 4- x26 + x27 + x28 >= 1 
xl7 4“ xl8 4" xl9 “H x20 4~ x21 4" x22 4" x23 x24 4" x25 “H x26 “b x27 -b x28 1 

x29 + x30 + x31 + x32 >= 1 

x21 + x22 4" x23 4" x24 4* x25 4* x26 + x27 + x28 4* x30 4" x32 >= 1 

x33 + x34 + x35 > = 1 

x31 "h x32 4" x34 “b x35 >=1 

x22 4" x26 4“ x27 H- x28 4* x32 4" x35 1 

END 

The solution to this linear integer prograraming problem has been obtained using LINDO 
software and is given in Table 5.3. The minimum production cost turns out to be $ 
2.061/piece. It is lower than that obtained using minimmn production cost/volume heuris- 
tic ($ 2.084/piece) and that obtained by adopting and the common practice of horizontal 
clubbing ($ 2.1223/piece). Thus there is significant improvement in production cost by 
using the suggested model which yields the minimmn cost. 

5.2.2 Example 2 

In this second example the geometry of stepped shaft is similar to that of the first example 
but the shaft dimensions are different. The shaft is shown in Fig. 5.2. 
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The mathematical model has been built on the basis of the same 15 machinable volumes 
and 35 candidate volumes. The minimum production cost for each candidate volume and 
other calculations are presented in Table 5.2. The minimization problem now gives 
minimize 

0.3859sl + 0.5821x2 + 0.8201x3 + 1.0275x4 + 0.2254x5 + 0.5149x6 + 0.7566x7 + 0.2725x8 + 
0.6492x9 + 0.9582x10 + 0.3848x11 + 0.4165x12 + 0.5545x13 + 0.6606x14 + 0.6495x15 + 
0.7952x16 + 0.1972x17 + 0.2443x18 + 0.3386x19 + 0.3976x20 + 0.2797x21 + 0.3268x22 + 
0.3268x23 + 0.4212x24 + 0.4801x25 + 0.3740x26 + 0.4683x27 + 0.5273x28 + 0.4880x29 + 
0.7183x30 + 0.5955x31 + 0.8716x32 + 0.1935x33 + 0.6009x34 + 0.7967x35 
subject to the similar 15 constraints as for example 1 because the shaft geometries are 
similar and the same candidate volumes have been selected in both cases. The Tuinirnmn 
production cost for the above problem turns out to be $ 2.5694/piece. It is important to 
note that the grouping of candidate volumes is different in this case even though the shafts 
are similar. Incidentally for this case the production cost using heuristic method and the 
common practice of horizontal clubbing turn out to be the same, i.e., $ 2.6353/piece. 

5.3 Remarks 

To sum up the mathematical model suggested in this chapter successfully gives the Tnim'mnm 
production cost for rough tm-ning of stepped shafts through optimal grouping of machinable 
volumes. The optimal grouping of machinable volumes achieved by mathematical model 
yields better minimum production cost than that obtained by using the common practice 
of clubbing these volumes horizontally or using heuristic method based on minimum pro- 
duction cost/volume. Similar stepped shafts may have different optimal grouping which 
depends on the size and location of machinable volumes. 



1 1 1 

2 1,2 

3 1,2,3 

4 1,2,3, 

5 5 5 

6 5,6 

7 5,6,7 

8 5,2 

9 5, 2,3,6 

10 5,2,3,4,6,7 

11 8 8 

12 8,9 

13 8,6 

14 8,6,3 

15 8,6,7,! 

16 8,6,3,^ 

17 10 10 

18 10,9 

19 10,9,7 

20 10,9,7,4 

21 10,12 

22 10,12,15 

23 10,9,12 

24 10,9,7,12 

25 10,9,7,4,12 

26 10,9,12,15 

27 10,9,7,12,15 

28 10,9,7,4,12,15 

29 11 11 

30 11,12 

31 11,14 

32 11,12 

33 13 13 

34 13,14 

35 13,14,15 


Table 5.1: Minimum production cost and production, cost/volume for compound volume. 



38 26 

40 26 

33 26 

38 26 

40 26 

38 33 

40 33 

38 33 

38 33 

40 33 

40 33 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

40 38 

38 28 

40 28 

38 28 

40 28 

28 17 


20.75 

20.75 

20.75 

20.75 

28.75 
28.75 
28.75 

49.5 

49.5 

49.5 


54.0 

54.0 
82.75 
103.5 
82.75 

103.5 
7.5 ■ 

61.5 
90.25 
111.0 

45.8 

66.5 

99.8 
102.55 

123.3 

120.5 
149.25 

170.0 
38.3 ■ 
38.3 

59.0 
59.0 


0.2177 

0.4172 

0.2537 

0.2797 

0.4796 

0.5246 

0.1580 

0.2153 

0.2458 

0.2678 

0.1986 

0.2206 

0.2559 

0.2589 

0.2809 

0.2779 

0.3084 

0.3305 

0.3914 

0.4072 

0.4408 

0.4652 
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SI. 

Seed 

Compound 

Do 

A 

dti 

L 

Qo 

Volume 

Cost /vol 

No. 

Vol. 

vohum^ 

mm 

mm 

mm 

mm 

$/ piece 

cu mm 

$ 10-ymm^ 

1 

1 

1 

26 

■B 

4.5 

50.0 

igH] 

15197 


2 


1,2 


II 

6.5 

50.0 

m 

23994 


3 


1,2,3 




50.0 

0.8201 

45357 

1.8081 

4 


1,2,3,4 



12.5 

50.0 


57923 

1.7790 

5 

5 

5 

30 

26 

2.0 

80.0 


14074 

1.6012 

6 


5,6 

38 

26 

4.0 

80.0 


48255 


7 


5,6,7 

42 

26 

8.0 

80.0 


68361 


8 


5,2 

30 

26 

2.0 

130.0 


22871 

1.1914 

9 


5,2,3,6 

38 

26 

4.0 

130.0 


78414 

0.8279 

10 


5,2,3, 4, 6, 7 

42 

26 

8.0 

130.0 

0.9582 

111087 

0.8536 

Bi 

8 

8 

38 

30 

4.0 

40.0 


17090 

2.2S;i8 

If 


8,9 

42 

30 

6.0 

40.0 


27143 

1.5345 

n 


8,6 

38 

30 

4.0 

120.0 


51271 

1.0816 

n 


8,6,3 

38 

30 

4.0 

170.0 


72633 

0.9095 

H 


8,6,7,9 

42 

30 

6.0 

120.0 


81431 

0.7977 

If 


8,6,3,4,7,9 

42 

30 

8.0 

170.0 


115359 

0.6893 

17 

10 

10 

42 

38 

2.0 

40.0 

0.1972 

10053 

1.9612 

.18 


10,9 

42 

38 

2.0 

80.0 

0.2443 

20106 

1.2151 

19 


10,9,7 

42 

38 

2.0 

160.0 

0.3386 

40212 

0.8421 

20 


10,9,7,4 

42 

38 

2.0 

210.0 

0.3976 

52779 

0.7533 

21 


10,12 

42 

38 

2.0 

110.0 

0.2797 

27646 

1.0117 

22 


10,12,15 

42 

38 

2.0 

150.0 

0.3268 

37699 

0.8670 

23 


10,9,12 

42 

38 

2.0 

150.0 

0.3268 

57805 

0.8670 

24 


10,9,7,12 

42 

38 

2.0 

230.0 

0.4212 

70371 

0.7286 

25 


10,9,7,4,12 

42 

38 

2.0 

280.0 

0.4801 

47752 

0.6823 

26 


10,9,12,15 

42 

38 

2.0 

190.0 

0.3740 

67858 

0.7832 

27 


10,9,7,12,15 

42 

38 

2.0 

270.0 

0.4683 

80424 

0.6901 

28 


10,9,7,4,12,15 

42 

38 

2.0 

320.0 

0.5273 

42223 

0.6557 

29 

11 

11 

38 

26 

4.0 

70.0 

0.4880 

59816 

1.1558 

30 


11,12 

42 

26 

8.0 

70.0 

0.7183 

66350 

1.2005 

31 


11,14 

38 

26 

4.0 

110.0 

0.5955 

93996 

0.8974 

32 


11,12,15,14 

42 

26 

8.0 

110.0 

0.8716 

8671 

0.9273 

33 

13 

13 

26 

20 

3.0 

40.0 

0.1935 

32798 

2.2316 

34 


13,14 

28 

20 

9.0 

40.0 

0.6009 

18778 

1.8321 

35 


13,14,15 

42 

20 

11.0 

40.0 

0.7967 

21314 

1.8592 


Table 5.2: Minimum production cost and production cost/volume for compound volumes. 
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SI. 

Compound 

IntcRor 

(’ornpound 

Min cost/ 

Min cost 

Compound 

Min cost 

No. 

volume 

prog. 

volume 

vol $ 10-5 

$/piece 

volume 

$/piece 

1 

1 

().;i3.5() 

8,6,3 

0.0969 

0.2797 

8,6,3 

0.2797 

2 

5, 2, 3, 6 

0.4330 

11,12,15,14 

1.2303 

0.4652 

12,15,10,19,7,4 

0.3305 

3 

8 

0.2177 

10,9,7,4 

1.9691 

0.2678 

11,14 

0.4408 

4 

10,9,7,4 

0.2678 

2,5 

2.4539 

0.3940 

' 2,5 

0.3940 

5 

11,12,15,14' 

0.4652 

13 

4.2459 

0.3417 

13 

0.3417 

6 

13 

0.3417 

1 

5.4796 

0.3356 

1 

0.3356 

■1 

total min 

2.061 





2.1223 

■ 

cost/piece 

math. 



heuristic 


common 

■ 


model 



method 


practice 


Table 5.3: Comparison of mimmuin production cost for the proposed method and common 
practice for stepped turning 


SI. 

Compound 

Integer 

Compound 

Min cost/ 

Min cost 

Compound 

Min cost 

No. 

volume 

prog. 

volume 

vol $ 10-5 

$/piece 

volume 

$/piece 

1 

8, 6, 3,4, 7, 9 

0.7952 

10,9.7,4,12,15 

0.6557 

0.5273 

8,6,3 

0.6606 

2 

11,14 

0.5955 

8,6,3 

0.7977 

0.6606 

12,15,10,9,7,4 

0.5273 

3 

10,12,15 

0.3268 

11,14 

0.8974 

0.5955 

11,14 

0.5955 

4 

2,5 

0.2725 

2,5 

1.1914 

0.2725 

2,5 

0.2725 

5 

13 

0.1935 

13 . 

2.2316 

0.1913 

13 

0.1913 

6 

1 

0.3859 

1 

2.5394 

0.3859 

1 

0.3859 


total min 

2.5694 



2.6353 


2.6353 


cost/piece 

math. 



hemistic 


common 



model 



method 


practice 


Table 5.4: Comparison of noinimum production cost for the proposed method and common 
practice for stepped turning 
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Chapter 6 

Conclusions and Scope for Future 
Work 

6.1 Conclusions 

The multipass turning problem has been solved using results of single pass turning problems. 
The three basic methods used to solve the single pass problem are: 

(i) Differential calculus with graphical search over nomographs of the constraints. 

(ii) Nonlinear programming (NLP) methods such as SUMT with DFP, GRG, requiring 
initial feasible points and methods which may be started from any arbitrary starting point, 
for example, SQP method, or which do not require any starting point such as some of the 
geometric programming (GP) methods. 

(iii) One dimensional search methods over tool Ufe and finding the maximum permissible 
feed imder the constraints for depth of cut. 

The last method has been selected due to simplicity in understanding and implementa- 
tion. The method requires simple calculations if the constraints axe monomials and it does 
not require any starting point. The parametric analysis single pass solution methodology 
has also been performed. It has been observed that optimal tool fife cutting speed and feed 
do not change with workpiece dimensions for turning operations provided the depth of cut 
and the cutting environment remain the same. The minimum production cost or time is 
obtained easily using simple calculations for changed workpiece dimensions. 

It has been a co mm on practice to remove a small depth in the finish pass and remaining 
depth of cut equally divided among all rough passes. An integer progr amm ing model has 
been developed on the basis of equal workpiece diameter. The model contains one fimsh 
pass and assumed mavimnm number of rough passes. It makes use of the production 
costs obtained for feasible series of depth of cut for finish as well as rough passes assuming 
some least count for depth of cut accuracy. The model has been solved using LINDO 
software. The optimal number of passes are also determined alongwith optimal depth of 
cut distribution. For each pass, depending upon the optimal depth of cut, optimal speed 
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integer programming model is logical selection of small depth of cut for finish pass alongwith 
optimal numbcm of passes which minimize the production cost per piece. 

It has Vx^en ol>served that the rough passes need not necessarily be of equal depth 
at minimxun pro<iuetion cost. This view has also been confirmed after solution of a few 
examples using intt'.gcr programming model. It has also been observed that sometimes 
LINDO software is not able to achieve or terminate at optimal results. Therefore, a second 
model based on. dynamic programming (DP) has been developed. Each pass represents a 
stage and each stage problem can be solved using any single pass solution methodology. 
Single pass problems are required to be repeatedly solved for a given depth of cut in each 
stage pass. The single stage/pass problem for the same depth of cut need not to be solved 
again and again if the single pass solution methodology suggested in the present work is 
used. The solution to new workpiece diameter and for the same depth of cut can easily be 
computed using parametric analysis. It has been observed that in cases where LINDO is 
not able to achieve oi.)timal results the dynamic programming approach clearly determines 
the optimal solution. 

The computation times of three nonlinear methods viz. SUMT with DFP, GRG, SQP 
and dynamic programming approach suggested in the present work have been compared 
for an example. The sequential quadratic programming approach treats all the passes 
simultaneously. The computer code-NCONF for this method has been obtained firom IMSL 
library and it may not be available for all PC users easily. This NLP method requires 
approximately the same order of computation time as is required for dynamic programming 
approach. The first two methods essentially require feasible starting point and considerably 
large computation time in comparison with DP approach alongwith single pass solution 
methodology. 

In the determination of optimal clubbing of volumes for rough turning of a stepped 
shaft, the common practice has been to club volumes horizontally only. In a similar work 
for milling, machinable volumes have been clubbed using a heuristic of minimum production 
cost per unit volume. The optimal clubbing of machinable volumes can be obtained using 
an integer programming model suggested in Chapter 5 for this purpose. First the feasible 
sets of compound volumes are formed and the corresponding minimiun production costs are 
obtained using DP method. Then these sets are chosen such that aU machinable volumes 
are covered and minimum production cost is obtained. 

To sum up there are following two approaches to solve multipass turning model : (i) 
solve all pas,ses simultaneously and (U) synthesize multipass solution using single pass solu- 
tions. The later one has been adopted in the present work because it decomposes complex 
multipass problems into simple single pass problems. The dynamic programming model 
has been used for this purpose. Out of various single pass solution methods viz. non-linear 
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programming techniques and graphical search over nomographs of constraints and chip- 
breaking region to find the optimal machining conditions, a simple one dimensional search 
over tool life has been selected to find the optimal machining conditions. At a given depth 
of cut, an optimal tool life can be located at which minimum production cost is obtained 
using the maximum permissible feed under the constraints. The optimal cutting tool life, 
speed and feed will remain the same for changed workpiece dimensions if the depth of cut 
and cutting environment is not changed. Dynamic programming alongwith single pass so- 
lution methodology has emerged as a good solution technique to obtained optimal cutting 
conditions for multipass turning with monomial constraints. 

6.2 Scope 

DP approach alongwith single pass solution methodology does not requires knowledge of 
any starting point. It also determines the optimal number of passes uniquely and is easy to 
implement. This may act as cutting conditions optimization module in computer assisted 
process planning of turning operations. The present work may be easily extended to obtain 
optimal machining conditions for plane and face milling operations where cutting tool life 
is s imil ar to extended Taylor’s tool life equation and constraints are monomial. It will also 
interesting to obtain the optimal machining conditions for continuous profile shafts using 
single pass methodology. 
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Appendix A 

SUMT with DFP Method 


Sequential Unconstrained Minimization Technique (SUMT) are penalty function methods 
which transform basic optimization problem into alternative opt imiz ation formulation such 
that the numerical solutions are obtained by solving a sequence of unconstrained minimiza- 
tion problems. Let the original optimization problem be stated as. 

Find decision vector X which minimizes 


/o(X) (A.1) 

subject to : 

Equality constraints 

/i(X) = 0 (A.2) 

where i = 1 . . . mj 

and inequality constraints 

/i(X) < 0 (A.3) 

where i = mi . . . m 

The above problem is converted to the following unconstrained optimization problem using 
interior penalty fimction method: 
minimize 


t—mi 1 t=m 1 

^(x, Ti) = /„(x) + g (A-4) 

If the minimization of (j) function is performed repeatedly for a sequence of values of penalty 
parameter Vj, the solution may be brought to converge to that of the original problem. The 
minimization of the (j) function has been performed using Davidon Fletcher Powell’s method 
discussed in the next section. • 
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A.l Davidon, Fletcher and Powell’s Method 


The iterative prorcHiiin^. of this method can be stated as follows. 

Step 1: Start with an initial point Xi and a n x n positive definite matrix Hi- Usually 
Hi is taken to b('. identity matrix I. Set iteration number ifc = 1. 

Step 2: Compute gradient of the function v/aj poiut and set 


Sk = -HkVA (A.5) 

Step 3: Find the optimal step length Aj^ in the direction of Sk and set 

Xk+i = Xk + A^Sk (A.6) 

Step 4: Test the new point for optimality. If Xk+i is optimal terminate the iterative 
process. Else, go to Step 5. 

Step 5: Update matrix Hk as 

Hk4-i = Hk + Mk + Nk (A.7) 


where 


and 


Mk = A 


* SkSk 

'sjQk 


(HkQk)(HkQk)^ 

QjHkQk 


Qk = V/(Xk+i) - V/(Xk) 


Step 6: Set the new iteration number — A: + 1, and go to Step 2. 


(A.8) 

(A.9) 

(A.10) 
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Appendix B 

Generalized Reduced Gradient 
Method 


Let the original problem be stated as 
minimize 



fo = MX) 

(B.1) 

subject to: 
equality constraints 

f(X) = 0 

(B.2) 

and bounds 

a<X<b 

(B.3) 


where X,a, b are N dimensional column vectors, f(X), is for any X under consideration, 
an m dimensional column vector and /o(X) is a real number. The fo and f are assumed 
continuously differentiable under the parallelotope defined by (B.3). The gradient of /o(X) 
is a row vector denoted by /g (X). Similarly, gradient of /j(X), i = 1 ... m is a row vector 
denoted by (X). The notation / (X) designates (m x n) matrix, the rows of which are 
= l...m. 

Vector X is partitioned into x and y, where y is m dimensional and x is n dimensional with 
n = N — m. Hereafter, y is called vector of basic or dependent variables and x is called 
vector of non-basic or independent variables. 

Notations (9/o/9x) and {dfo/dy) are called row vectors of fo with respect to x and y. 
Similarly, (m x N) matrix /’(X) is partitioned into (m x n) matrix (df/dx) and (m x m) 
matrix {df/dy). 

Let a feasible point be partitioned as to satisfy the following nondegeneracy assumption. 
There exists a partition of X into x and y such that 

a; < y] < bj'ij (B-4) 
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(9f/9y°) stands for {df/dy) computed at x° and y° . 
The GRG algorithm can be summarized as follows: 


Step 1: Compute — hi®, the direction of move for independent variable x, using the following 
substei>s. 


Step 1 . 1 : Compute the reduced gradient 



Step 1.2: Compute the projected reduced gradient (p°) given by \fj, 


(B.5) 


.,0 

Pj 


0 if aj° = aj and —g^ < 0 
0 if = bj and —gj > 0 
g^ otherwise 


Step 1.3: Compute h°, the modified projected reduced gradient, i.e., the (opposite) gra- 
dient direction of move for the independent variable x. This direction may simply be 
h.° = p°. However, it has been modified by making all the components such that — %) 
or — bj) is too small and equal to zero. It has been checked that after these modifica- 
tion scalar product (g°) x (h°) is positive. 


Step 2: Compute a first value of positive number 9i.- 

Step 3: Compute x° — 0ih°, and project it onto parallelotope aj < Xj < bj.'ij, to obtain 
X i.e. set 


Vi, 


X = 


aj i/ if x° - 0ih° < aj 

bj i/ifx°-0ih°>6y 

if x° — 0ih° otherwise 


Step 4 : Compute a feasible corresponding to 0i, i.e., to solve with respect to y the 
system of m equations in m unknowns: 

f(x\y) = 0 (B.6) 


This is done by some iterative method described after step 5. 

Step 4 .I: If no speedy convergence is achieved then decrease 9i = 0-56i and go to step 3, 
with same h®. 

Step 4.2: Otherwise, let y^ be the solution obtained for eqns. (B.6) and the corre- 
sponding point in N dimensional space. 
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Step ) - /o(^°)5 decrease $i as above and go to step 3 with same h°. 

Step 0(;h<'rwis<^, at the end of step 4, we have some feasible which satisfies 

/o(X^) < /o(X°) (B.7) 

Step 5: At; this st,a.g('. cnther set X“ := X^ and begin a new iteration, or try to improve the 
last obtained valu<‘ lor (9i. Ixi doing so step 3 is reached for any improved value of 9i, with 
same h® and (Wi'ut.ually terminate with some X^ satisfying eqns. (B.6) and then begin new 
iteration, X® :=' X'. 


B.l Iterative Method and Change from Basic to Non-basic variables 


The solution to eqns. (B.6) means a solution 


|f(x^y)| = e 


(B.8) 


where e is a given positive number (tolerance), and | . | is some norm in m dimensional 
space. Let , y^ and so on upto y* be successive iterates obtained in appljdng step 4 with 
in the same iteration of GRG. The term ‘no speedy convergence occurs’ implies that eqns. 
(B.6) does not hold for any t < to- Consequently , the iterations in step 4 are stopped as 
soon as eqns. (B.6) is satisfied for some y*. 

Set y^ = y*, to obtain X^ in step 4-^- A. pseudo-Newton method has been used in our code 

. . . . 

(B.9) 


= y‘ 




If = (x^,y*) € P and = (x^,y*'''^) B P, then intersection of segment X* — X*"^^ 

with boundary of P is calculated. For some index r yj. = or b^. The basic variable y^ 
is changed with non-basic (independent) variable Xj for with some index s such that the 


following expression is maximized Vj 
min ( I arj (&; - xj) , | (xj - aj)) 

where | a,.;, is element of row and colmnn of the matrix (^) (^). 
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Appendix C 


Sequential Quadratic Programming 
Method 


NCONF subroutine is based on subroutine NLPQL, a FORTRAN code developed by Schit- 
tkowaski [1986]. It uses a sequential quadratic progr amming method to solve the general 
nonlinear programming problem. Let the original optimization problem be stated as 
minimize 

/o(X) (C.l) 

subject to : 
equality constraints 

/i(X) = 0 (C.2) 

where i = 1 . . . mi 

and inequality constraints 

/i(X) > 0 (C.3) 

where z = mi . . . m . It should be noted that the code NCONF has been made to accept 
greater than or equal to form of the constraints. Therefore all the constraints in multipass 
turning problem which are in less than or equal to form have been modified accordingly 
and 

a < X < b (C.4) 

where all the functions involved in the problem have been assumed to be continuously differ- 
entiable. The method based on iterative formulation and solution of quadratic progr ammi ng 
subproblems , obtains subproblems by using a quadratic approximation of Lagrangian and 
by linearizing the constraints. That is, 
minimize 

1 /2d^Hfcd -h ( vf(X;t))^d (C-5) 

subject to : 

V d + /i(Xt) = 0 (C.6) 
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